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ABSTRACT 


Design of multivariable systems is considered and design pro- 
cedures are formulated in the light of the most recent work on model 
matching. The word model matching is used exclusively here to mean 
matching the input-output behavior of two systems . In the frequency 
domain the term is used to indicate the comparison of two transfer 
matrices containing transfer functions as elements. The design of 
multivariable systems is particularly complicated by the fact that the 
transfer matrix is a ratio of two polynomial matrices. The use of 
state variable feedback not only affects the pole positions, but also 
alters the positions of most of the zeros of individual transfer 
functions due to interaction. Because of these complexities, non- 
interaction has been one of the main criteria of design in the past. 

This study concentrates on design methods where non-interaction 
is not used as a criteria. Non-interacting design is obtained as a 
special case of the more general interacting design. Two design methods 
are considered. In the first method the design is based solely upon 
the specification of generalized error coefficients for each individual 
transfer function of the overall system transfer matrix. The main 
disadvantage of such a design is that it does not take into considera- 
tion the transient response of the system, as the concept of error co- 
efficients is based solely upon steady state behavior. Thus, the 
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transient response of the system may not be satisfactory. However, if 
the transient response is satisfactory, the main advantages of the 
first method is that it is simple and gives a realizable solution to 
the problem. 

The second design method is called the pole fixing method 
because in it all the system poles are fixed at preassigned positions. 
In addition, the zeros of terms either above or below the diagonal 
can be partially fixed via steady state error coefficients. This is a 
considerable improvement over the conventional decoupled design where 
it might not even be possible to decouple the system. Even if the 
system were decouplable, some of the system poles could be uncontrol- 
lable. However, designing the system using the pole fixing method 
requires that certain sufficiency conditions be met. In many cases 
these sufficiency conditions may be satisfied by a name change in the 
output vector. 

An example is worked to demonstrate the use of both design 
methods. The special case of triangular decoupling and minimum con- 
straints are discussed. 
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INTRODUCTION 

The problem of designing a linear multivariable system has long 
occupied the minds of control engineers, especially recently in view of 
its vast application in engineering and related fields. Most of the 
earlier work is centered on designing noninteracting systems. Efforts 
started as early as 1950 and continued through most of the 1960 ’s 
[Boksenbom and Hood, 1949; Kavanagh, 1956, 1957, 1958; Horowitz, 1960; 
Morgan, 1964; Rekasius , 1965], Falb and Wolowich [1967] first formu- 
lated and proved necessary and sufficient conditions for determining 
whether or not linear state variable feedback (l.s.v.f.) can decouple 
a multivariable system. They also gave a design procedure for placing 
the poles of the decoupled system at desired places. Then Gilbert 
[1969] gave a physical interpretation to the results of Falb and Wolo- 
wich and explained the notion of feedback invariance. Gilbert's work 
provided the complete answer to the problem of decoupling a multi- 
variable system using l.s.v.f. 

But the problem of design of multivariable systems was far 
from solved, except for the special case of decoupling. Luenerger 
[1567] devised a transformation matrix which puts the state equations 
of multivariable system into a canonical form, similar to phase vari- 
able form for the scalar case of single input-single output systems. 
Wolowich and Falb [1969] made use of the canonical transformation 
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devised by Luenberger [1967] and came up with startling results on the 
structure of multivariable systems. 

While a lot of attention was devoted to using dynamic compensa- 
tion for multivariable systems [Wolowich, 1973; Moore and Silverman, 
1971; Morse and Wonham, 1970; Pearson and Ding, 1969; Silvinsky, 1969] 
studies continued utilizing l.s.v.f. for design of interacting control 
systems [Ferg, 1971; Heinz, 1968; Anderson and Luenberger, 1967]. The 
design of interacting multivariable systems is given the name Model 
Matching, In view of model matching some pertinent questions arise. 

Does there exist a ,set of feedback and gain matrices {F,G} such 
that when applied to a given system the input-output behavior of the 
system characterized by the frequency domain transfer matrix T(s) 
matches a pre-specif ied input-output behavior. Wolowich [1972] 
answered this question in great detail. He utilized the structure 
theorem [Wolowich and Falb , 1969] and for the first time developed a 
complete set of feedback invariants for multivariable systems. He also 
gave necessary and sufficient conditions for Exact Model Matching, and 
outlined the procedure to find the pair {F,G),if the necessary and 
sufficient conditions were met. 

However, the main design question of what to do if the 
necessary and sufficient conditions are not met remains unanswered. 
Unfortunately, this is the most important problem a designer is likely 
to face. The necessary and sufficient conditions for model matching 
do not give any hint as to what can be expected from the system, or 
what changes need to be made in the model transfer matrix to ensure 
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that the necessary and sufficient conditions are satisfied. The designer 
is thus left with no guidelines with which to proceed. 

In this study a completely new and basic approach is taken to 
the design of multivariable systems. The approach provided at least 
a partial answer to the above questions. The approach taken in this 
study is based upon generalized error coefficients. The design speci- 
fications for multivariable systems could either be expressed directly 
in terms of generalized error coefficients or could be specified as 
transfer functions. In the latter case the transfer function is 
represented in terms of its Maclaurin series so that the coefficients 
of the successive powers of s have a direct relationship with the 
generalized error coefficients. 

A design based solely on error coefficients may not be satis- 
factory for the two reasons. First, only a specified number of error 
coefficients can be incorporated into the design. Secondly, the 
error coefficients describe the system input-output behavior only 
after sufficient time has elapsed for system transients to decay to 
Insignificant amplitudes. Since the pole positions are not known 
until after the design is complete, decaying transients are not ensured 
in advance . 

A second design technique called pole fixing design method takes 
care of the above uncertainty by pre-specifying all the pole positions 
in advance. The main feature of this method is that in addition to 
arbitrarily fixing the poles, it maintains enough freedom to specify 
error coefficients of one side of diagonal terms in the transfer matrix, 
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thus in turn, indirectly specifying system zeros . This method has the 
following distinct advantages over the conventional decoupled design, 

1. All the off-diagonal terms are not forced to zero, 

2. There are no uncontrollable poles, and 

3. Those systems which cannot be decoupled could very well be 
designed using pole fixing method.- 

A disadvantage of the pole fixing method is the requirement that 
the system satisfy certain additional sufficiency conditions. However, 
these sufficiency conditions are variant under the name change of 
inputs and outputs and are also dependent upon the output matrix of 
the system. Thus, in many cases ,the sufficiency conditions could be 
satisfied by simple changes in the system output matrix. 

This study has been organized as follows. In Chapter 2 nota- 
tions are introduced and a brief treatment of background material is 
provided. This includes a discussion of the feedback invariants and 
structure of the multivariable systems. The representation of a trans- 
fer matrix as a ratio of two polynomial matrices is used to derive the 
central result of this study in Chapter 3, The relation between the 
error coefficients and the components of polynomial matrices represent- 
ing the transfer matrix is obtained. This relation produces a series 
of equations which can be evaluateji in sequence. 

Additional notations are given in Chapter 4 as a means of 
denoting the submatrices of a component matrix. An identity is intro- 
duced to partition the multiplication of two matrices into multiplica- 
tion of its submatrices. The design freedom is illustrated with the 
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help of Theorem 4.1. The proof of . the theorem results in a 
step-by-step design procedure for the design of multivariable systems 
utilizing the generalized error-coefficients as the design criteria. 

An example is used to illustrate the design steps. 

Chapter 5 develops. the pole fixing design method. The short- 
comings of the error coefficient design method are overcome by the pole 
fixing method. Once again the design constraints are introduced via 
Theorem 5.1. A step-by-step design procedure is given for a quick and 
easy reference and the example of Chapter 4 is reworked using pole 
fixing method. 

This study is concluded with Chapter 6 where the results of 
this study are summarized and further research is suggested. 



CHAPTER 2 


STRUCTURE OF MULTIVARIABLE SYSTEMS 

2.1 Introduction and Organization of the Chapter 
This chapter provides a review of existing work on the structure 
and exact model matching of multivariable system. Section 2,2 describes 
the system representation in state variable form followed by some defini 
tions and notations. Section 2.3 describes a transformation [Luenberger 
1967] which produces a phase variable representation for multivariable 
systems. The above transformation is used in Section 2,4 where the 
structure theorem of Wolowich and Falb [1969] is explained and proved. 
The advantage of such a structure is that it separates the system input 
output behavior transfer matrix in two parts, one which is invariant 
under state variable feedback and, second, which is almost completely 
dependent upon state variable feedback. The above property is used to 
derive necessary and sufficient conditions for exact model matching 
[Wolowich, 1972], as briefly discussed in Section 2.5. 

The past work on the decoupling of multivariable systems has 
intentionally not been discussed because almost all the information 
about decoupling is provided by the structure theorem. Hence, in 
Section 2.6 the problem of decoupling is dxscussed as a special case 
of the problem of exact model matching. It is shown hew all the 
relevant information about decoupling can be obtained from the use of 
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structure theorem and from the necessary and sufficient conditions for 
exact model matching of the multivariable systems. 

Finally, this chapter is concluded with a discussion of the 
advantages as well as hazards of attempting to use the above mentioned 
necessary and sufficient conditions for model matching. ^ 

2.2 System Representation, Definitions and Notations 
Throughout the study it is assumed that the linear multivariable 
time invariant open loop system is represented in state variable form 
by following well-known equations. - . 

X = AX + BU 

Y = CX (2.1) 

where 

X G system state vector 

Y c ^ system output vector 

U e system input vector 

The matrices A, B, and C are constant system matrices of appropriate 
dimensions. The following state variable feedback is considered 

U = FX + GR (2.2) 

where F is mxn and G is mxm constant matrices and R is new reference 
control input vector. With the state variable feedback as specified in 
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Eq. (2.2), closed loop state variable representation of the system is 
given by 

♦ 

X = (A + BF)X + BGR 

Y = CX . (2.3) 

While the system is described completely in the time domain, 
much of the work contained here has been transformed into frequency 
domain. The following definitions and notations prove useful. 

Definition 2.1 Transfer Matrix : A transfer matrix is any 

matrix whose elements are transfer functions representing 
the input-output behavior of a single input - single output 
system in the frequency domain. 

The transfer matrix associated with .the closed loop system of Eq. (2,3) 
is denoted by T(s) and is given by 

T(s) = C(si - A - BF)’^BG (2.4) 

Definition 2.2 Proper and Strictly Proper Transfer Matrices : 

A transfer matrix is called proper (strictly proper) if for 
every transfer function element, the degree of the numerator 
polynomial is less than or equal to (strictly less than) the 
degree of the corresponding denominator polynomial. 

The transfer matrix of the closed loop system of Eq. (2.3) given by 
Eq. C2.4) is a strictly proper transfer matrix because there is no 
direct feed forward from inputs to outputs. 
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Definition 2,3 Polynomial Matrix : A polynomial matrix is 

any matrix whose elements are polynomials of degree greater 
than or equal to zero, 

A polynomial matrix Is usually denoted by capital letter with explicit 
S dependence unless otherwise mentioned. The transfer matrix of 
Eq, (2.4) can be described as a product of two polynomial matrices as 
T(s) = R(s)P(s) Here P(s) ^ represents the inverse of the polynomial 

matrix P(s). 

Unless otherwise mentioned, the i^^ row and i^^ column of a 
constant matrix is denoted by a subscript i and superscript i, respec- 
tively; viz., denotes i^^ row of B matrix and b^ deonotes j^h 
column of B matrix. The 1^^ row of a non constant matrix is denoted 
by a subscript on capital letter, viz., P^(s) denotes 1^^ row of P(s). 

The following notations are equivalent for any constant or 
polynomial matrices. 


B = [b^^] 


P(s) [p.^ (s)] 

-^J 


(2.5) 


where b , denotes the element at the cross section of i^“ row and 
column of B matrix and p^^ (s) denotes polynomial at the cross section 
of row and i^^ column of polynomial matrix P(s). The determinants 
of. a constant matrix B and a non constant matrix P(s) are denoted as 


follows 
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determinant B = {Bj 

determinant P(s) - |p(s)l (2.6) 

2»3 Luenberger Canonical Transformation 
Various canonical transformations are available for multivari- 
able systems. The one given by Luenberger [1967] is particularly useful 
because the matrices representing the transformed system can be parti- 
tioned into submatrices which have a phase variable form. The property 
of the phase variable representation of a system is that for the scalar 
case of single input-single output system, all the relevant information 
about the input output behavior of the system is contained in the last 
row of system matrices. The system transfer-function can be written 
by inspection of system matrices. 

Even though the above does not hold in the multivarible case, 
i.e., the transfer matrix cannot be written by inspection, the trans- 
formation is useful because it separates the system transfer matrix in 
an invariant part and another part which arbitrarily depends upon state 
variable feedback except for its form. This is similar to the scalar 

i 

case where it is well known that zeros of transfer function are invariant 
and the poles can be arbitrarily placed by state variable feedback. 

Only the order of system cannot be increased. 

To f:‘nd the transformation matrix and the transformed system, 
consider the system described in Eq. (2.1). It is assumed that the 
input matrix B has rank m (otherwise, the m inputs are not indepen- 
dent) and the output matrix C has rank m (otherwise, the m outputs are 



11 


not indopsTident^ . Xt is nlso sssuinod. tliat the system described in 
(2.1) denoted by {A,B,C} is controllable, i.e., the nxnm controllability 
matrix 

[B: AB: ... A^”^B] (2.7) 


has rank n. 

That is to say the above matrix has at least n independent 
vectors; however, in general there may be more than one set of n inde- 
pendent vectors. 

The following algorithm finds a transformation matrix denoted 
by Q which has the properties as described above. 

Step 1 ; Find the first n independent columns of the controllability 
matrix given in Eq. (2.7). Let the B matrix be represented by 


B ^ [b^ ... b“] 

where b^ represents i^^ column of B matrix. 

Step 2 : Rearrange the n independent columns found in Step 1 in the 

following manner 


[ b 9 Ab 9 b^Oy u ^ bj 


where as indicated A^b^ is the column of nxm matrix A^B. Notice 
that the above rearrangement groups ' together o^ independent vectors 
which are related to i^^ column b^ of the input matrix B. Thus, 
informally speaking, the n^^ order multivariable system has been 
divided into m subsystems with the i subsystem having order . 
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Notice that 


^ 1 for i = 1 . . . m. 


and 


ID 

I <J. = n' . 

i=l 


Step 3 : Define the integers and the Ixn. row vectors as follows: 


K 

dj^ = I for K = 1,2, ..., (2.9) 

i=l 


Z = d th row of the inverse of the matrix 

Jx JX 


in Eq. (2.8) 


Step 4 : Form the nxn transformation matrix Q as follows 


Z^A 


Z^A 


01-1 


Q = 


Z^A 


V 


02-1 


Z A 
m 


Om-l 


(2.10) 
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A .1 

If the transformation X = QX is used, where Q is the nonsingular 


A 


matrix as given above and X is 

the new state vector. 

then the new system 

equations are given by 



II 

A A A 

AX + BU 


Y = 

A A 

cx 

(2.11) 

where 



H 

O' 

II 

A = 


B = QB s 

B = 


t = cq“^ 

C = CQ 

(2.12) 

and the state variable feedback 

is given by 


A A A 

U = FX + GR 


(2.13) 

where 



FQ = F and 

A. 

G = G 

(2.14) 


As mentioned above , the main advantage of such a transformation 
is that the matrices A and B take a special fonn as shown below » 



(2.15) 
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where for i=j is a x matrix in phase variable form as shown 
below. 


11 


0 10 

0 0 1 

• 

« 

0 0 0 

-k it it 


0 0 

0 0 

0 1 

it it 


(2.16) 


where **s denote non zero terms in general. 

For 1 i , A. , is a. X a. matrix with all but its a. th row iden- 
ij 1 J i 

A 

tically equal to zero. The matrix B in Eq, (2.8) takes the following 
special form. 


A 

B = 


0 


di,l "di,2 


do,l do, 2 




d-i ,m 


d o ^m 


m* 


(2.17) 


where are as defined in Eq. (2.9). Thus, only d„ th rows of B are 

K. IS. 


non zero. 
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It is worth noticing at this point that the transformation 
matrix Q obtained in Eq. (2.10) is not the only transformation matrix 
which places A and B in the form given in Equations (2.15) and (2.17) 
respectively. In general there may be other arrangements for picking 
n independent columns out of the controllability matrix of Eq. (2.7) 

A A. 

which lead to a transformation matrix so that A and B have the above 
mentioned forms. That is to say, in general the o’s (roughly speaking 
is the order of i^^ subsystem) are not unique, but in the absence 
of any other general method for finding the transformation matrix Q, 
the algorithm described in this section is used throughout this study. 


2.4 Structure Theorem of Wolowich and Falb 


To be able to better understand the structure theorem it is 
necessary to give a few definitions introduced by Wolowich and Falb [1969] 
in connection with the structure theorem. 

Definition 2.4 is defined as the mxm diagonal 

matrix as given below. 



(2.18) 


Thus the i^^ diagonal entry of the matrix 3^ is 



16 


Definition 2.5 S (s ) t S(s) is defined as the following nxm 
matrix of single terra monic polynomials. 


S(s) - 


,n-l 


='^2-1 


,cfm-l 


S^(s) 0 0 ... 0 

0 S^(s) 0 ... 0 


0 ... S^(s) 
m 


(2.19) 


where 


S.(S) = 


,^i-l 


is a X 1 column vector. 


Definition 2.6 A: A is defined as raxn matrix consisting of ,the 


m ordered d th rows of A = Thus, A is given as fellows: 

K 




A = 




^ d ^,2 


^ 6-2 ^ d 2 » 2 






do ja 




( 2 . 20 ) 
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Definition 2.7 B is defined as the tnxm matrix consisting 

of the m ordered d^^ the rows of % = Thus, B is given 

as follows: 


K 


ij 


B = 


di,l X,2 


d2,l d2,2 


d-, ,m 


d ^ , in 


b j bj A ... bj 


( 2 . 21 ) 


Equivalently if a matrix E is defined as follows 


E = 


0 0 
0 0 


. 0 

, 1 

. 0 

« * 

• « 

• • 

0 0 


* « 


C2.22) 


Then clearly 


B = EB 


(2.23) 


In view of above definitions and the form of transformed 
matrices A and B as given in Equations (2.3.5) and (2.17), respectively. 
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^ A A A. 

notice that BF has all but the m ordered th rows zero. Thus, A + BF 

Iv 

A. 

has exactly the same form as A. Let 


$ = [d>. .] = A + BF = 
ij 


^11 ^12 **• ^Im 


ml m2 mm 


(2.24) 


Then for i=i , <5^. is a. x cr. matrix in phase variable form as 
ij 1 3 

shown in Eq. (2.16) and for i j , is a x matrix with all but 
last row identically zeros. Thus, i has exactly the same form as A as 
given in Eq. (2.15). 

^ .iw A 

Hence, as in Eq. (2.20), define $=A+BFasa mxn matrix 

consisting of the m ordered d th row of the matrix ^ as defined above 

K 

2 

in Eq. (2.24) and partition it into m row vectors as follows: 



HA' 

$ 

11 

^12 

$ ' 
Im 

? = A + B F = 

% 

21 

?22 ••• 

22 


ml 

¥ 

m2 

$ 

mm 


(2.25) 


Then, is 1 x row vector obtained by partitioning the d^ th row 

of into m row vectors. Thus 


^^il ^i2 *** ^im^ ^d^,2 '** 


(2.26) 
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The structure theorem is now derived which essentially separates 
the system input-output transfer matrix into invariant and dependent 
parts with respect to state variable feedback. As mentioned earlier in 
the Definition (2.1), for the multivariable system of Eq. (2.1) with 
feedback given in Eq. (2.2), the input-output transfer matrix is given 
by Eq. (2.4). Substitution for A, B, and G from the set of Equations 
(2.12) in Eq. (2.4) gives 

T(S) •= "^Q(sl - Q”^^Q - Q'^BF)”^ Q^^BG 
= C(sl - ^ - BFQ~^)BG 
= C(sl - A - Bf)”^ §G 

^ A A 

Substitution for A + BF from Eq, (2.24) in above gives 

T(s) = C(sl - ^)~^ BG (2.27) 

where 

= F q“^ (2.28) 

It is claimed by Falb and Wolowich [1969] that the system trans- 
fer matrix of Eq, (2,28) can be written as follows, 

T(s) = C S(s) [(Sg)”^ {S® - ? S(s)] 


(2.29) 
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or, equivalently 


C(sl BG = C S(s) [S^ - $S(s)]“^ BG (2.30) 

where S^, S(s), are as defined in Eqs, (2.18), (2.19), and (2.25), 
respectively. In their proof of the claim, B is not as general as given 
by Eq. (2,17). The proof given below is more clear and is consistent 
with the notion introduced in this section. 

Post-mult iply both sides of Eq . (2.30) by G Then, to prove 
Eq. (2.30), it suffices to prove the following; 

v> 

(si - $)’^ B = S(s)[S^ - is(s)]"^B 

After a trivial manipulation, it can be shown that it suffices to prove 
the following: 


(si - $) SCs) = g B ^ IS^ - ‘^S(s)] 

After substituting for % from Eq. (2.23), it now becomes necessary to 
prove only the following: 


(si - t) S(s) = E[s‘’ - 4S(s)] (2.31) 

Expansion of the left hand side of the above equation results in 
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(sI-<I>)S(s) = 


-sMii 

12 

_$ 1 
Im 

21 

• 

« 

sI-$ 

22 

— $ 

2m 

« 

ml 

m 2 

. . . sl-^ 

nsi-fij^)Sj^(s) 

-♦ 1282 ( 3 ) 


(s) 

(sI-4>22)S2(s) 


(s) 

. -*m2®2(®> 


S^(s) 

0 



0 

S2(s) 



0 


m 




(sI-®^)S (s) 
mm m 


( 2 . 32 ) 


Now examine the (i,j)th block matrix in the above equation for i“j 


CsM^^)S^(s) = 


s ~1 
0 s 


0 

-1 


0 0 0 

* * * 


s 

* S-(|) 


0 

0 

-1 

di,di 


.''i-2 


s 


ai-1 



0 

. S. (s) 
1-3 ^ 



22 


zero. 


The above is a cr^ x 1 vector with all but the last element non 
For i j 




Again, this is a x 1 column vector. 

After substituting for these block matrices in Eq. (2,32), the 
following is readily obtained. 
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(sI-$)S(s) - 


^llS^Cs) 


■^ 12^2 




"'* 21 ®!^®^ 


■ * 22 ® 2 <®) 








- $ s (s) 
nnn m 


' s ^ ~ (s) -4>-^S-(s) 










s^ - $ S (s) 
mm m 


where E is as defined in Eq, (2.22), 
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The above Eq. (2.33) can be written as follows: 


(sI-^)S(s) = E 


S^l 

0 

0 

* * » 

0 

s ^ 

0 

* • • 

• 

♦ 

♦ 

• 

9 

0 



11 

12 

• • « 

Im 

21 

$ 

22 

« • « 

% 

2m 

* 

• 

^~ml 

m2 

• • • 

% 

mm 




Sl(s> 


S^(s) 


ta 


} 


and substituting for S° and S(s) from Eqs. (2.18) and (2.19), respec- 
tively, one gets 

(sl“^) = - '^S(s)} 

The above is exactly the result claimed in Eq. (2.31). Thus, the claim 

A 

of Ealb and Wolowich [1969] is demonstrated to be true for B of 
Eg. (2.17), which is more general than the one used by them. 

It is thus established that system input output transfer matrxx 

T(s) can be written as follows: 


T(s) = C S(s) [(BG)'^(S° - *S(s))r^ 
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After substituting froin Eq. {1.15) for ^ one obtains, from above 

^ ^ n * ^ ^ •^1 

T(s) = C S(s) [(BG) is - (A + BF) S(s)}] 

that is, 

r(s) =■ R(s) P(s)"^ (2.34) 

where 

R(s) = C S(s) (2.35) 

is a mxm polynomial matrlK independent of state variable feedback and 
thus forms the invariant part of transfer matrix for any state variable 
feedback of the type shown in Eq. (2.2), and 

PCs) = [(BG)"^ (S^ - (A + BE) S(s)}] (2.36) 

is also a mxm polynomial matrix. Since P(s) is a function of both F 
arid G, P(s) takes into account all the effects of the state variable 
feedback of Eq. (2.2). 

It is worth noticing at this point that P(s) not only accounts 
for all the feedback and gain, but 'has a unique form for a given set of 
sigmas. Recall the sigmas are fixed (even though not unique). Also, 
notice that the coefficient of the highest order term in each (ij)th 
polynomial is given by the corresponding (i,j)th term in the matrix 



(BG) Thus, by choosing G properly, coefficients of highest order 

tenn in each polynomial of P(s) can be arbitrarily manipulated. 

The coefficients of all other lower degree terms in the poly- 
nomials of P(s) can be arbitrarily placed by proper choice of feedback 
matrix F because B is a nonsingular matrix. 

Thus, the conclusion is that, except for its form, the poly- 
nomial matrix P(s) is almost completely arbitrary. The only exception 
is that the matrix coefficients consisting of the highest degree poly- 
nomial in each column of PCs) must be a nonsingular matrix because that 
matrix is (BG) ^ which is a nonsingular matrix. The above property of 
a polynomial matrix is referred to as being column proper by Wolowich 
[1972]. 

2.5 The Necessary and Sufficient Conditions for 
Exact Model Matching 

The term model matching is more easily explained in the scalar 
case of single input-single output systems. Given a scalar system 
(plant) and a scalar model if by some means (here by state variable 
feedback) , the plant input-output transfer function can be made equal 
to the model transfer function, then the system (plant with feedback) 
is said to be matched with the model because their input output behav- 
ior is identical. In the scalar case it is well known that 

1) The zeros of plant transfer function are feedback invariants 

2) The order of the closed loop plant can not be increased by 


state variable feedback. 
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3) If the plant is controllable, then the poles of the plant 
transfer function can be arbitrarily placed by using 
state variable feedback. 

4) The closed loop gain can be adjusted at desired level. 

Knowing the above feedback Invariants and dependents, the necessary and 
sufficient conditions for exact model matching are easily obtained for 
the above scalar case and are well known. Namely, the system can be 

matched to the desired closed loop transfer function, that is, the 
model, by state variable feedback alone if, and only if, 

1) The model has zeros at the same place as that of the plant. 

2) The pole zero excess of the model is the same as that of 
the plant . 

3) If the model has zeros other than plant zeros, they must 
cancel with the model poles. 

In the scalar case the knowledge of the above is necessary and 
sufficient conditions for model matching lead directly to a design pro- 
cedure. Even though the terra exact model matching for multivariable 
system is thought of as just an extension of the scalar case, it is 
actually more than this, i.e., instead of matching two transfer functions, 
one has to match two transfer matrices. Several additional definitions 
are needed to extend the discussion to the multivariable case. 

Definition 2.8 Unlmodular Matrix ; A polynomial matrix is called 

unimodular if the determinant of a matrix is a non zero polynom- 
ial of- degree zero, i.e., is a scalar. 
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Definition 2.9 Division of a. Polynomial Matrix : A polynomial 

matrix RCs) is said to be right [left] divisible by another 
matrix HCs) if R(s) H(s) ^ [H(s) ^R(s)] is also a polynomial 
matrix. (Consequently , H(s) is called right (left) divisor 
of R(s) 0 

Definition 2.10 Relatively Prime Polynomials ; Two or more 
polynomials (taken together) are called relatively prime, if 
there does not exist any polynomial of degree greater than zero 
which divides all the polynomials simultaneously. 

Definition 2.11 Relatively Prime Polynomial Matrices : Two 

polynomial matrices are said to be relatively right [left] 
prime, if there does not exist any non-unimodular matrix which 
right [left] divides both the matrices, (Notice that any uni- 
modular matrix divides any polynomial matrix and, hence, while 
making a test for relative primeness only non-unimodular 
matrices may be considered.,) 

The necessary and sufficient conditions for exact model matching 
are given in the following two theorems which are due to Wolowich [1971, 
1972] . 

Theorem 2.1, A (mxm) rational strictly proper transfer matrix 
can always be factored as the product. 

T(s) = R(s) P“^Cs) 


(2.37) 
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where RCs) and P(s) are relatively right prime polynomial 
matrices, respectively, in the Laplace operators, with P(s) non 
singular and column proper, 

A formal constructive proof of this theorem is given in [Wolowich, 1971] 
and is not repeated here. 

Theorem 2.2, Consider the system given in Eq. C2.1) with 
{A,B} controllable and B of full rank m < n, and let the 
the model transfer function be given by 

T (s) = R Cs) P"^(s) (2.38) 

m m m 

where andP^C®) are relatively right prime polynomial 

matrices. There exists a linear state variable feedback 
(l.s.r.v.) pair {F,G}, with G non singular, which satisfies the 
relationship : 

T(S) = C(sl - A - BF)“^BG = R(s) P"^(s) 

= R^(s) P"^(s) = T (s) (2.39) 

m m m 

if, and only if, for some non singular polynomial matrix, H(s) 
the following three conditions hold ; 

1) R(s) = R (s)H(s), l.e. , R(s) is left divisible by R (s). 

m m 

2) The m ordered a. of P (»)H(s) are identical to those of P(s), 

i m ' ' 

3) Pjjj(s)H(s) is column proper. 
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The proof of sufficiency is constructive, and simply involves equating 

P (s)H(s) to P(s), once an appropriate H(s) has been found, i.e,, from 
m 

Eq. (2.39) and substituting for P(s) from Eq, (2.36) 

P„(s)H(s) = PCs) = CBG)"^ IS° - (A + BF) S{s)] (2.40) 

m 

Since P^(s)H(s) is column proper, the mxm matrix consisting of the 
coefficient of highest degree s terms in each column of P^(s)H(s) is 
non singular. Let this matrix be D then 

D = (BG) 
or 

G = (2.41) 

The corresponding F can then be determined by first premultiplying 
Eq. (2.40) by (BG) , and then substracting from both sides, i.e., 

(BG) P (s)H(s) - S° = - (A + BF)S(s) (2.42) 

Since the m ordered o’s of P^(s)H(s) are the same (condition 2) 
as that of P(s)i in the left hand side of above Eq. (2.42), a complete 
cancellation of occurs, and the left hand side could be written as 
-M S(s) where M is some constant matrix. Hence 

-M S(s) = ~(A + BF)S(s) 

The above is satisfied if the following holds 



31 


or 


A + BF = M 
F = b“^(M-A) 


or 

F = FQ = {b"^ CM-A)} Q 


(2.43) 


Thus, the set {F,G} Is uniquely determined for a given H(s), which satis- 
fies the conditions of Theorem 2.2. 

To establish necessity, consider Eq. (2.39); since P(s) is non- 
singular 

R(s) = R (s)P"^(s) P(s) 
m m 

= R^(s) P*(s) PCs)/ |f^(s)| C2.44) 

where P*(s) is adjoint P (s) and ]p (s) 1 = det . P (s) . 
m m ' m m 

Since R (s) and ? (s) are relatively right prime, it follows that 
m m 

[MacDuffee, 1956] there exists two polynomial matrices and M^Cs), 

such that the following relation is satisfied. 

M, (s)R Cs) + M-(s) P (s) = I (2.45) 

I m 2 m m 

where I is the mxm identity matrix, 
m 

Postmult ip lying both sides of Eq, (2,45) by P(s) , one 


obtains 
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M, (s)R (s)P*(s)P(s) + M_(s)p (s)F*Cs) P(s) = P*(s) P(s) (2.46) 
1 m d z m m m 


But from Eq. (2.44), |p (s) I divides R (s)P*(s)p(s) because R(s) is a 

' m mm 

polynomial matrix; also because 


P„(s) p*(s) = |p„(s)| I 
in m 'mm 

|P^(s)| divides both left side members of Eq. (2.46), which implies it 
must also divide the right side of Eq. (2.46). 

It can thus be concluded 


. P*(s)P(s) 

P J^(s) P(s) = 

m 


is a nonsingular polynomical matrix, since it has been established 


P (s) I divides P*(s)P(s). But from Eq. (2.44) 
m ’ m 


R(s) = R (s) [P“^(s)P(s)] 
m m 


= R (s) H(s) 
m 


(2.47) 


where H(s) = P^^(s)P(s) is a polynomial matrix. It is thus proved that 

R(s) is left divisible by R (s) which is condition 1 of the theorem. 

m 

Also since 


H{s) = P^'-(s) P(s) 
m 


(2.48) 


P (s)K(s) = P (s) 
m 


(2.49) 
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Equation (2,49) directly implies conditions 2 and 3 of the theorem and, 
hence, the theorem is proved. 

In summary, finding a pair {F,G} for exact model matching con-- 
sists of the following steps: 

(a) finding the transformed representation of the system. 

(b) factoring the desired transfer matrix (model) T^(s) as a 

product of two relatively prime polynomial matrices 

R (s) P‘^(s). 
m m 

(c) determining the appropriate H(s), if one exists. 

(d) determining the pair {F,G} by comparing P^(s)H(s) with P(s), 

There are algorithms available to perform the first two of these steps. 
If the system is invertible, i.e., if R(s) has rank m, then H(s) , if it 
exists, is uniquely given by Eq. (2.48). 

However, note that just employing the above equation does not 
always result in a proper H(s) which satisfies all three conditions of 
Theorem 2.2. Thus, in general, there is no algorithm to find a suit- 
able H(s), but if the system is Invertible, then all one can find is 

f 

whether there exists or not an H(s) which satisfies all the three condi- 
tions of Theorem 2.2. 

It is clear from the above discussion that the necessary and 
sufficient conditions for exact model matching, as given in Theorem 2.2, 
car: be successfully employed to test whether a given system could be 
matched to a given model (even though the test is cumbersome). However, 
unlike the scalar case, the necessary and sufficient conditions fail to 
predict what can be matched to the system. In terms of the design no 



hint is given as to what changes are required to satisfy the necessary 
and sufficient conditions if a design being evaluated fails to meet 
these conditions* The only case where these questions are completely 
answered, is the case of a decoupled model, which is discussed next. 


2.6 The Problem of Decoupling 

Necessary and sufficient conditions for decoupling a multivari- 
able system were first given by Falb and Wolowich [1967] and were inter- 
preted by Gilbert [1969] ; The above work has been examined in detail by 
others as well [Slivinsky, 1969; Ferg, 1971; Agrawal, 1972], and, hence, 
no attempt is made to repeat it. Rather, it is assumed that the multi- 
variable system in consideration is decouplable, i.e., necessary and 
sufficient conditions for decoupling are satisfied. Let the following 
decoupled transfer matrix be realizable and nonsingular. 


T^(s) = 


dii(s) 


n (s) 
mm 

d (s) 
mm 


(2.50) 


where is relatively prime to d^^(s) for i = 1 . . . m. 

The above assumptions imply that there exists a pair {F,G} which 
realizes Eq. (2.50) and the structure theorem of Section 2.4 implies that 


T.(s) can be written as follows: 
a 
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where R(s) and P(s) are as given in Eqs. (2.35) and C2.36), respectively. 
Post multiply both sides of Eq. (2.51) by P (s) to obtain 


T^(s)P(s) = R(s) (2.52) 

Now consider the ith row of both sides of Eq. (2.52) which is given by 


i.e. , 


n. . (s) 

0 . . . 0, /g\ > 0 • • • 0 


n (s) 


P(s) = R^(s) 


(2.53) 


where as mentioned previously in Sec. 2.2, a subscript of i on a matrix 
denotes ith row of the matrix. 

Let us assume P^(s) and r^(s) to be the greatest common divisor 
of the elements in P^(s) and R^(s) respectively, and let P|(s) and R|(s) 
be the prime polynomial vectors left after taking out the greatest com- 
mon divisor polynomials p^(s) and r^(s) . Then the following relation is 
obtained , 

P^Cs) = P^(s) P^(s) 

R^(s) = r^(s) R^(s) (2.54) 

Substituting for P^(s) and R^(s) from above in Eq. (2.53) and multiplying 
both sides by d^^(s) the following result is obtained. 

{n,,(s) p,(s)} P!(s) = {d,,(s)r.(s)} R! (s) (2.55) 

ii i 1 ii 1 1 

Notice that {n^^(s) p^(s)} is a polynomial multiplying a prime polynomial 
vector P^(s) and similarly, ^ polynomial multiplying 



36 


prime polynomial vector R^Cs)« only way the above is possible is 

that 

Pj^(s) = k R[(s) 
which from Eq. (2.55) implies 


HiiCs) P^Cs) k = r^(s) 


(2.56) 


where k is a scalar nonzero constant. From above Eq. (2.56), one 
obtains the following; 


k n^^(s) r^(s) 

p.(s) (2.57) 

xi 1 

Now since n^^(s) is relatively prime to d^^Cs) , there must exist 
a polynomial fCs) of a degree greater than or equal to zero, such that 

r^Cs) = k fCs) 

Pi(s) = diiCs) f(s) C2.58> 

From Eq.(2.57) and Eq. (2.58), it is clear that all the zeros 
of subsystem [zeros of n^^(s)]are contained in zeros of r^(s) . In 

addition, there may be other zeros of r^(s) [zeros of f(s)l which cancel 
with zeros of p^(s). To show that zeros of r^(s) need not cancel with 
zeros of p.(s), i.e., all the zeros of r,(s) do appear as zeros of i 
subsystem, unless otherwise intended, it only need be shown that there 
exists a pair {F',G} such that 
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r^(s) k 

d.,(s) f’Cs) 


(2.59) 


where p^(s) = f^Cs) is relatively prime to r^(s) . 

Let us choose f’(s) such that it is relatively prime to f(s), 
but the highest degree term in both polynomials is exactly the same. 
Thus, the coefficient of highest order term in p’Cs) is the same as that 
in p^(s), and the degree of p^(s) is equal to the degree of p^(s). Now 
form the vector polynomial 


P^(s) = p|(s) P^(s) 


where P!^(s) is defined in Eq. C2.54), namely 

P^Cs) = p^(s) P|(s) 


(2.60) 


(2.54) 


Equations (2.60) and (2.54), together with the fact that p^(s) and 
p.(s) have the same degree and the same coefficient of the highest degree 
term, imply that PV(s) has the same form as P^(s). Then, if P(s) is 
column proper, the P'*(s) is also column proper, where 


P”(s) 


P”(S) - 
L J 


That is to say, if there exists a pair {F,G} which realizes P(s), then 
there exists a pair {F* ,G} which realizes P”(s). 
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It is thus proved that Eq. (2,59) holds, i.e,, zeros of r^(s) 
need not cancel with the zeros of p V(s) . Since it has already been 
shown that zeros of subsystem are contained in zeros of » it 

now implies that zeros of r^(s) are the zeros of subsystem unless 
they are intentionally cancelled with the poles of subsystem. 

Since R(s) is Invariant, ^nd, hence, r^(s), is invariant. 

It can thus be concluded that zeros of r^(s) are invariant zeros of 

subsystem, if decoupling is to be maintained. Also, all the poles 
of i^^ subsystem can be arbitrarily placed because p^(s) can be made an 
arbitrary polynomial except for the coefficient of highest degree term 
without changing the form of P(s) and without affecting its property of 
column properness. 

The actual order of i^^ subsystem is easily obtained because 
the degree of P^^(s) equals the sum of the degree of r^(s) and the pole 
zero excess of the i^^ subsystem, both of which are knovm. 

The only other information that remains to be obtained concerns 
the so-called uncontrollable poles. It is well-known that if the sum of 
the orders of all m subsystem does not equal n, the order of the overall 
system, then there are some poles of the system which are unaccounted 
for. Gilbert [1969] and Silvinsky [1969] called them uncontrollable 
poles because these unaccounted poles cannot be controlled by state 
variable feedback, if exact decoupling is to be maintained. 

In practice, seldom is the model or the feedback so accurate as 
to be able to force the system to be exactly decoupled. Thus, these 
uncontrollable poles appear together with the other subsystem poles 
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and only in case of exact decoupled system do they cancel with system 
zeros at exactly the same position. Thus, If any of these uncontrollable 
poles are in the right half s plane, then the system is an unstable sys- 
tem, There is nothing one can do about these poles if decoupling is to 
be preserved because these poles are feedback invariant. So to ensure 
a stable decoupled system, one must know the positions of these so-called 
uncontrollable poles. It is shown next that the position of these poles 
is known from the knowledge of R(s), the invariant part of the system 
tranfer matrix, without even actually decoupling the system. . Consider 
Eq. (2.5l) snd take the determinant of both sides 

= |r(s)p”^(s)| = |r(s) | |p"^(s)| 

= |r(s)| / |p(s) 1 
or 

m , . 

d ~ (sT “ |R(s)1 / |P(s)l C2.61) 

1=1 

Now, from Eq. C2.54) ^j^Cs) = r^(s) R’ (s) , and hence. 


R(s) = 


r^(s) R|(s) 


r (s) R’ (s) 
m m' 
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and 


where 


|R(s)| 



r^(s) 1 r'Cs)1 


R’(s) 


R{(s 

R'(s) 

m 


(2,62) 


Substituting the above expression for jR(s) | 


xn ii.q. , 




gets 


IT 


1=1 


d..(s) 


TT 1R'(s)1 

1=1 


1p(s)| 


Substituting for 


nii(s) 


from Eq, (2.57), one gets 


or 


m 



i=i 


r (s) 

1 p^(s) 


r^(s) |r*(s)| 
i=l 

1p(s)1 


K. 


m 

TT Pi(s) 

i=l 


P(s) 


(2.63) 


where Kj[ is some scalar constant. 
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Since the left side numerator is a nonzero constant in Eq. 

(2»63), the right hand side numerator polynomial jR'(s)| must divide its 
denominator polynomial jp(s) j , Thus, zeros of [r* ( s) [ must be zeros of 
|p(s)|. But zeros of |r’ ( s) j are Invariant under state variable feed- 
back (because R(.s) is invariant and, hence, R'Cs) is invariant) and 
zeros of jp(s)| are the overall system poles. That is to say, zeros of 
}r^(s)| form the invariant poles of overall. system. These Invariant poles 
are the so-called uncontrollable poles of the overall decoupled system. 

Thus, these uncontrollable poles are actual poles of. the over- 
all system and are cancelled by zeros in the same position if exact 
decoupling can be obtained. However, as mentioned earlier, in practice 
seldom is the model or feedback and gain so accurate as to produce 
exactly the decoupled system. Thus, rarely is the pole zero cancella- 
tion exact. That is to say, if any of these uncontrollable poles are 
in the right hand s plane, then the system is bound to be unstable if 
decoupled. 

Since R(s) is completely known |r' (s) j is known and hence 
positions of uncontrollable pole can be found in advance without com- 
pleting the actual design of the decoupled system. If any of these 
poles are in right hand s plane, one might as well give up the idea of 
decoupling the system. However, it should be remembered that the so- 
called uncontrollable poles are uncontrollable only as long as de- 
coupling is required. Otherwise, if the overall system is controllable, 
all the n poles of the overall system can be arbitrarily placed and 
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thus, a completely stable system could be obtained if decoupling is not 
a criteria, 

2.7 Summary and Conclusions 

The main objective of this chapter has been to provide a compre*- 
henslve review of past and recent work on multivariable systems. Follow- 
ing the comments on system notation and a few definitions, Luenberger's 
[1967] canonical transformation is described to represent multivarible 
systems in phase valuable form. Next, it is shown how Wolowich and Falb 
[1969] used this transformation to find a structure for multivariable 
systems. The main advantage of such a structure is that it separates 
the system input output transfer matrix into state feedback invariant 
and state feedback dependent polynomial matrices. This property is used 
to find necessary and sufficient conditions for exact model matching . 
Lastly, the special case of decoupling is discussed and it is shown how 
the structure of multivariable systems in conjunction with necessary ' and 
sufficient conditions for model matching can be used to find all realiz- 
ability information about decoupled models. 

However, except for the decoupled models, it has been shown that 

the necessary and sufficient conditions for exact model matching fail to 

produce any realizability criteria. Thus, the necessary and sufficient 

conditions are useful only for testing whether a given model can be 

realized from the system. Even the test irj cumbersome and difficult 

computationally, because it requires that the model transfer matrix be 

put into the form T (s) = R (s) P ^(s) where R (s) and P (s) are 
^ m m m mm 
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relatively right prime polynomial matrices. Thus, if the test fails, 
then the designer has to start all over again to find some other model 
for which the test conditions might be satisfied. However, he has no 
guidance as to what changes to make and no assurance that the new model 
might be any more realizable than the previous one. 

In the next chapter an important result is obtained which is 
subsequently utilized in the following chapter to give completely new 
and step-by-step design procedures for the design of multivariable 
systems. 



CHAPTER 3 


THE DESIGN EQUATION 

3«1 Introduction and Organization of the Chapter 

This and the following chapters form the main contribution of 
this study. The results of this chapter form the backbone for the 
design methods presented in the next chapter. 

In Section 3,2 the concept of generalized error coefficients is 
extended to the multivariable case. The transfer matrix T(s) and the 
two polynomial matrices R(s) and P(s), which are state variable feed- 
back invariant and dependent respectively, are- represented in terms of 
their component matrices. The notation for the derivatives of these 
matrices is also introduced. 

In Section 3.3 a result is obtained which relates the deriva- 
tives of the polynomial matrix R(s) with the derivatives of the transfer 
matrix T(s) and the polynomial matrix P(s). The next section, 3.4, 
utilizes this result to give a relation among the component matrices of 
R(s)j P(s), and T(s). This relation is the main result of the chapter. 
Lastly, a theorem is proved which is vital to the application of the 
main result as described above to the design of multivariable systems. 
The theorem leads directly into the next chapter. 
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3.2 The Generalized Error Coefficients 
for Multivariable Systems 

For the scalar case of the single input-single output system, 
the generalized error coefficients are directly obtained from the 
Maclaurin series expansion of the transfer function of the system 
[Truxal, 1955] as given below. Let t(s) be the transfer function of a 
single input-single output system and y(s) and u(s) be the Laplace trans- 
forms of the output and input, respectively, then 


y(s) 

u(s) 


= tCs) 



(3.1) 


In the classical and original definitions 

= position error constant 

= velocity error constant 

k = acceleration error constant 
a 


The above error coefficients are measures of the steady state errors if 
the input is a unit step, unit ramp, unit parabolic function, and so on, 
respectively. 

The primary disadvantage of the classical definitions rests in 
the limited amount of information available from the specification of 
error constants , since only one constant is significant. The general- 
ized error coefficients represent an attempt to circumvent this diffi- 
culty by defining all error constants in terms of the low frequency 
behavior of y(s)/r(s). With the Maclaurin series expansion of t(s) as 
shown in Eq. (3.1), the error constants are defined in terms of succes- 
sive coefficients in the series. 



46 




Equation (3.1) could be written as follows: 

y(s) = u(s) ^ s 'i(s) ^ ^(s) - (3.2) 

p V a 

If the inverse transform of series of Eq. (3.2) is taken term by term 
and impulses at t=o are neglected, y(t) is given by 

k ^ 

y(t) = ^ u(t) - — u’(t) - — u"(t) . . . (3,3) 

p V a 

where dashes indicate derivatives with respect to time. 

When the transfer function t(s) is replaced by its Maclaurin 
series in Eq. (3.1), the transient terms of y(t) are discarded. Thus, 
Eq. (3.3) is a valid description of the output y(t) only after suffi- 
cient time has elapsed to allow those terms in y(t) which are generated 
by poles of t(s) to decay to insignificant amplitudes, furthermore, the 
validity of Eq. (3,3) evidently depends upon the rapidity of convergence 
of the series of Eq. (3.1) [Savo, 1953]. 

Thus , only if care is exercised in the use of generalized error 
coefficients can it be said that the generalized error coefficients 
describe the relation between output and the reference input. 

Equation (3.1) can be written in a more convenient manner as 

follows: 

2 

t(s) = k^ “ k^s - k^s - . . . 


(3.4) 
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Clearly, 




and so on. 

In the case of multivariable system the transform of the output 
vector Y(s) is related to the transform of input vector U(s) as follows: 


Y(s) = T(s) U(s) (3.5) 

where T(s) = [t^^ (s)] is an mxm transfer matrix with the (ij)th entry 
being the transfer function t^^(s). 

Since every entry in the transfer matrix T(s) is a transfer func- 
tion, a Maclaurin series expansion of the form of Eq. (2.4) could be 
written for each (i,j)th entry as follows; 


ti^(s) = k.j, 



(3.6) 


With the expansion of each term of T(s) as in Eq. (3.6), T(s) can be 
written as follows 


T(s) = - K^s 



K s 

r 


r-1 


(3.7) 


where 


K = [k., ] 

r lir 


is an mxm constant matrix. Thus, is called the rth component of 
transfer matrix T(s). 
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From Eq. (2.34), T(s) = R(s)P(s) where R(s) and P(s) are 
polynomial matrices defined in Eq.. (2.35) and (2,36), respectively. 
Thus, R(s) and P(s) can also be written in the same form as that of 
Eq. (3.7) as follows: 


R(s) = + R 2 S + ^2® + • . • 


P(s) = Pj^ + P^s + P^s + . . . 


(3.8) 


where R^ and P^ are ith component matrices of R(s) and P(s), respectively 
The main difference between Eq. (3.7) and (3.8) is that while 
Eq. (3.7) is an infinite series, Eq. (3.8) is a finite series. Let: 


” {T(s) } be denoted by 


{R(s)} be denoted by R^ 


(3.9) 


{P(s)}be denoted by P^ 


P ^ ^ 

r+1 “ r ! 


r-1 r! 
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In the next section a relation is obtained between the rth 
derivative of R(s) and the derivatives of T(s) and PCs). 


3.3 The Derivatives of the Feedback 
Invariant Matrix R(s) 

Consider Eq, (2.34) of the last chapter 


T(s) = R(s) P"^(s) 


(2.34) 


Post multiply both sides of Eq. (2.34) by P(s) to obtain 


R(s) = T(s)P(s) 


(3.10) 


In view of Eq. (3.10) and the notation of Eq, (3.9) the following result 
is proved next 


r’' = I t 


(3.11) 


1=0 


for 


r > o 


where 


r ! 


(r-i)! i! 


Proof: The result is proved by induction. Clearly, for the result is 

true for r=o because for r=o , Eq. (3.11) becomes 


R = CT P=TP 
o 
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i.e. 


^ 

ds° ds° ds° 


(s) 


or 

RCs) = T(s) P(s) 
which is true from Eq. (3.10). 

Let Eq. (3.11) be true for the r^^ derivative, then expanding 
Eq. (3.11) for r^^ derivative. 


= r_ T^P° + r + . . . 

L* 

o 1 


+ r, p" + r^ + . . . 

^i ^i+1 


+ r^ + r 

t> - Vi 

r-1 r 


(3.12) 


Differentiate the above Eq. (3.12) with respect to s to obtain the 
following; 


+ r P^ + r P^ + 

C Vj C., o- 

o oil 


^1 1+1 


+ r_ ^r-i-lpr+a ^ ^ ^ . + r + r . 

^i+1 ^r-1 r-1 


1 r o r+1 

+ r,. TP + r T P 

O Vi 
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Combining the terms with similar derivatives in T(s) and PCs), one 
obtains 


+ Cr^. + 

o o 1 


But , 


+ . . . (rc ^ ) T-" P^^^ 

^ ^i+1 


+ . . . (r^ +r^ ) P*" + T° 

r-1 r r 


- 1 _ 

^ ^o 

o 


^r+1 

r 


( 3 . 13 ) 


and 


+ r 


r! 


rr 


i!(r-i)! (i+1) ! (r-i-1) ! 


r ! 


(i+1) ! (r-i) ! 


{i+1 + r-i} 


(r+1) ! 


(i+1) J (r-i) ! 


r+1 


‘^i+i 



After substituting from the above identities in Eq 


(3 . 13) , the 


following result is obtained: 


o 1 i+1 


+ p’^ + T°P’^+^ 

r r+1 


or 


, - r+1 , - 

^r+l ^ j r+1 

. ^ 1 
r=o 


^r+l-ipi 


(3,14) 


Comparing Eq. (3.11) ith Eq. (3.14), it is evident that Eq. (3.14) is 
the same as Eq. (3.11) with r replaced by r+1. 

Thus, the result of Eq, (3.11), namely 


R 


r 


i=o 


r T 
C. 

1 


r-i 


(3,11) 


has been established by induction for any integer r ^ o. 

In the next section the central result which links component 
matrices of T(s) with component matrices of R(s) and P(s) is established. 

3.4 Relation Among Component Matrices 
In this section a relation is obtained among the various component 
matrices of T(s) in Eq, (3,7) and the component matrices of R(s) and P(s) 
in Eq. (3.8). This relation is examined to produce a design procedure 
and a theorem concerning necessary and sufficient conditions for the 
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After evaluating the derivative of Eqs . (3.7) and (3.8) at 
s=o (steady state) and after using the notation of Eq. (3.9), one gets 


s=o 


•^1 


for r = o 
for r > 1 


R 


s=o 


iRl 

r!R 


r+1 


for r=o 
for r > 1 


s=o 


1 

r!P 


r+1 


for r = o 
for r > 1 


(3.15) 


Now, evaluate Eq, (3.11) at s=o and substitute from Eq. (3.15) to 
obtain 


= I rc T 
s=o i=o i 


r-x 


s=o 


s=o 


Now, evaluate the above equation for r=o to obtain 


Rl = 


(3.16) 


and evaluate for r ^ 1 to obtain the following result 


r ! 


o_r 


s=o 


= rc TP 
r 


s=o 


r-1 

+ I r r 

i=o ^i 


:-x 


s=o 


for r > 1 


s=o 
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After substituting for the derivatives from Eq. (3.15), one 


obtains 


r-1 


r!R 


r+1 




i+1 


x=o 


1 .e. 


T“1 

r!R = r!K,P - T r! K .P. 
r+1 1 r+1 r+l-x x+1 

i=o 


for r > 1 


Cancelling the common non zero coefficient r! from both sides of the 
above equation and taking the summation term on the left, the result 


IS 


r-1 


K,P _ = R _ + I K ,P_, for r > 1 

r+l-i x+1 — 


1 r+1 r+1 


x=o 


The above can be x^ritten as follox«7s 


r-1 

.1 == R .1 + I K . P. 

1 r+1 r+1 r+l-x x+1 

1^-IU 


+ K P 
r+1 1 


(3.17) 


Break Eq, (3.17) in two parts, one for r=l and the other for r ^ 2, and 
write together with Eq. (3.16) as follows: 
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for r=o 


for r=l + K 2 ^1 


r~l 

for r > 2 K,P . = R ,- + y K . P, 

— 1 r+1 r+1 r+l-i i+1 

1=1 


r+1 1 


(3.18) 


The set of Eqs . (3,18) could be written as one single equation 
for all value of r > o as follows: 


K, P ^ = Q + K 

1 r+1 r r+1 1 


(3,19) 


where 




\ “ ^ 1^1 


R, 


for r = 0 
for r = 1 


r-1 

+ I K _ . P_, for r > 2 
r+1 r+l-i i+1 — 


i=l 


Evaluating Eq. (3.19) for r = 0, 1, 2, 3 . . .in sequencej the following 
sequence of equations is obtained: 


''i \ 

= R2 + 


K, P, = R, + K-P. + K,P, 
i J i I c J i 


(3.20) 
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From the sequence of Eqs. (3.2), it is evident that 
. . . can be calculated in sequence for any given ^2* ^3 * ’ * 
respectively, if in nonsingular. Thus, nonsingularity of is a 
sufficient condition for unique solution to the components of P(s). The 
following theorem is presented here to prove the necessity of being 
nonsingular for any useful design. 

Theorem 3.1: K^, P^, and all must be nonsingular for 

any stable design for which transfer matrix TCs) is 
nonsingular. 

Note that the reverse is not true; i.e. , the system could be unstable 
and still have nonsingular K^, and . However, the theorem tells 
that if any of P^, or is singular, then either the design is 
unstable or the transfer matrix obtained is singular. 


Proof of Theorem 3.1: From Eq. (2.36) P(s) can be calculated for any 

given {F,G}. Also, from Eq . (3.8) 


Then 


Let 


= ?(s) 


s=o 


= !P(s) 


s=o 


|p(s)| = (s + x^)(s + . is + X^) 


Then 
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1pCs)1 


s=o 




X 

n 


or 






Now, if the system design is a stable design, then none of the 
eigenvalues Is zero and, hence, 4 - 0; i.e., is nonsingular. Note 

once again that the above does not imply that X^s are all negative or 
have real part negative. All it implies is that if is singular, then 
one of the eigenvalues is zero, resulting in an unstable design. 

Now, from Eq. (2.34) 

T(s) = R(s) P ^(s) 
or 

P(s) = t“^(s) R(s) 

since T(s) is assumed to be nonsingular. 

From the above, one easily obtains 


PCs) 


Let 


= t“^Cs) 

R(s) 

s=o 


s=o 

T"^(s) 




s=o 



s=o 


then from above, it follows that 
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?! 


But the above equation and the fact that must be nonsingular implies 
must also be nonsingular. However, evaluating Eq. (2.34) at s=o 
gives 


T(s) 

= R(s) 

p'^(s) 


s=o 

s-o 


i.e. , 

The above is valid because is nonsingular; i.e., P^^ exists. 
Now, since both matrices on the right are nonsingular, it implies 
must also be nonsingular. Hence, the theorem is proved. 


3.5 Summary 

The Information developed in this chapter is new and forms the 
basis for the design of the multivariable systems as discussed in detail 
in the next chapter. The main result obtained in this chapter is the 
identity of Eq. (3.19) from which a step-by-step design procedure can 
be developed. The sequential nature of design procedures to be developed 
is indicated by Eq, (3.20), where it is seen that each component of P(s) 
can be determined from the knowledge of the corresponding component of 
T(s) and previous components of P(s) . 

In the last section a theorem (Theorem 3.1) is presented which 
established usefulness of Eq. (3.19). If existence of inverse were 
not necessary for a stable design v;ith the nonsingular transfer matrix, 
then one would have put too much of a restriction on K, by saying that 
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must be nonsingular. Ideally, one would like to see all elements in 
to be unity because (ij)th element of is nothing but the value of 

y^(s)/u (s) 

^ s=o 

i.e. , steady state response to a unit step input. But, by Theorem 3.1, 
it is seen that if is singular, then either at least one of the sys- 
tem pole lies at origin or the transfer matrix is singular (the output 
responses are not independent) both of which are undesirable situations. 
Thus, Theorem 3.1 gives significant information about what steady state 
response can be achieved for multivariable systems. 



CHAPTER 4 


A DESIGN METHOD FOR MULTIVARIABLE SYSTEMS 

4.1 Introduction and Organization of the Chapter 
In this chapter the results of the previous chapter are utilized 
to develop a design method and to indicate constraints on design 
requirements. First, additional notation is introduced in Sec. 4.2 
to compactly represent columns, rows and elements of the component 
matrices. Next in Sec. 4.3 a' theorem concerning the design constraints 
is proved. The proof of the theorem is constructive. A step-by-step 
design procedure is outlined in Sec. 4.4. The design procedure is 
illustrated with the help of a simple example in Sec. 4.5. In the last 
section the results of the chapter are summarized and the advantages and 
drawbacks of the design method are discussed. It is pointed out how 
some of the unwanted features of this design method could be overcome. 
This is done in the next chapter . 

4.2 Notation and Identity 

Unless specifically mentioned, the following notation is used 
throughout . 

1. A polynomial matrix P(s) is represented in the following 
equivalent manners . 

P(s) = Ip^^(s)] = ^ 3 ®^ + ... + **' 
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so that p . • (s) denotes (ij)th element of the matrix P(s) 

J 

and 


P 


r+1 


1_ 

rl 


~ PCs) 
ds^ 


s=o 


is a component matrix of PCs) for r = 0, 1, 2 ... 

A kxA submatrix of a component matrix P^ is denoted by speci- 
fying all its elements inside a kxil constant matrix followed 
by a subscript r as shown below; 


^ij ^ij+1 
P- . 1 - 

XtIj 






^1+kZ 


^i+k j+£ 


r 


Also Cp. .) denotes Cij)th element of P 

1 J f r . 


The ith column of a component matrix P^ is denoted by the 
following equivalent notations: 



Thus, a group of columns of P^ can be denoted in the fo] lowing 
equivalent ways: 
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[Pr Pr 


i+1 


j 1 r i i +1 j 1 

p^] = [p p ••• p ] 


Clearly, the above is a matrix with j-i+1 columns in it. 
Accordingly, can be denoted in terms of its columns as 
follows : 


T> r 1 2 

P = [p p 

r *^r ^r 


m, .12 
P^i = IP P 


m, 
- P ] 


4. In an earlier section, a subscript is used to denote the row 

of a matrix. Here, however, subscripts are used for components 
of polynomial matrices. Thus, to avoid confusion and to be 
consistent with an earlier notation, the component matrix or 
the submatrix of a component matrix is first denoted by a 
single capital letter. The rows are then denoted by the usual 
notation, i.e., by a subscript. For example, to denote the 
rows of a component matrix, p^ , let W = . Then denotes 

the ith row of the matrix P^ . Thus, the notation introduced 
here is compatible with the earlier notation where a sub- 
script i denotes ith row and a superscript 1 denotes ith column 
of a constant matrix. 

Next , an identity is established which is useful in proving the 
theorem of the next section. The identity essentially partitions the 
the multiplication of two matrices in such a way that the known part is 
separated from the unknown part. 
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Identity A.l: Let VW = Z where V, W, and Z are constant 

matrices of dimension mxp, pxq and mxq, respectively Cci, 

p, and q are integers). Let a. th, a^th ... a th denote the 

X z p 

p ordered columns of V and the same p ordered rows of W. 

If j m, then the product Z = VW can be broken into two 
matrices as shown below. 


Z = [v“l v“2 



+ v“j+2 ... 

"^“j+1 




L 


_ 


(A.l) 


That is to say, Z is the sum of two matrices, the first of which 
is made from j ordered columns of V and corresponding j ordered rows of 
W. The second matrix is the multiplication of the remaining (p-j) 
ordered columns of V and the corresponding (p-j) ordered rows of W. 

The above in effect is a way of partitioning the multiplication of two 
matrices into two parts. 

Proof of the Identity: Clearly 


z = VW = [v^ v^ . . . v^J 


w. 


w. 


w 

p 


( 4 . 2 ) 



where, as stated above, v denotes ith column of V and w, 

3. 

denotes ith row of W, The above Eq. C^.2) can be rewritten 
as follows: 

V W- + V W- + V w_ + ... = Z 

12 3 p 

Thus, Z is the sum of p matrices each of dimension (mxq) as 
shown above. By reordering, the above can be written in the following 
equivalent manner. 

P 

Z = ^ v^w = v^'^ w + v*^^ w + ... w 

^ “1 «2 


From Eq. (4.3), by combining the p reordered columns of v and the p 
ordered rows of W, one gets: 


Z = [v“l v“2 . 


OIt^ t 
V Pj 


w. 


ai 


w. 




w 




(4.4) 


From Eq. (4.4) it is clear that if the p columns of V are reordered 
and the p rows of W are reordered accordingly, then the result of the 
multiplication is unchanged. 


Now, partition Eq. (4.3) as follows: 
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“j+1 , °^j+2 

+ V + V + ... V ‘ 


'aj+i "'ctj+2 




(4.5) 


In the above equation combine the first j matrices in one batch and 
the remaining (p-j) matrices in the second batch to obtain the required 
identity of Eq. (4.1). The identity is thus established. This identity 
is used in the proof of Theorem 4.1. 


4.3 Design Method I 

The design described in this section is based mainly upon the 
result derived in the last chapter, especially as expressed in Eq. 

(3.19). It is seen from the form of Eq. (3.20) that ... can 

all be determined in sequence from the specification of *»•> 

respectively. It is also seen that the generalized error coefficients 
describe relations between the output and the reference input. However, 
it is unfortunate that in the general case the must satisfy 

additional constraints, and thus the K^’s are correspondingly restricted. 
The constraints on P(s) are briefly described in Section 2.4 and 2.5. 

To begin the development, Eq. (2.36) is now examined in view of 
necessary and sufficient conditions for eract model matching as given 
in Section 2.5. Recall that in Section 2.5, P(s) is given as 
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PCs) = CBG)“^{S^ - (A + BF) S(s)} 

= CBG)'*^ - (BG)'"^ (A + BF) S (s) (2.36) 

Irrespective of what B and G are, if (BG) exists, then any 

polynomial in the ith colixmn of P(s) has a maximum degree of This 

-- _X 

is true because any polynomial in the ith column of (BG) (A + BF) S(s) 

has a maximum degree of cr^-1, and any polynomial in the ith column of 

~1 o ^ — 1 

(BG) S has a maximum degree of o^. However, (BG) must exist for 

realizability, since the condition that P(s) be column proper for model 

~ -1 

matching comes from the existence of (BG) 

Thus, for realizability the degree of the highest order poly- 
nomial in each column of P(s) must remain unchanged as this has been 
shown to be one of the necessary and sufficient condition for model 
matching in Sec. 2.5, Also P(s) must be column proper for a nonsingular 
transfer function (Sec. 2.5). By definition P(s) is column proper if 

the matrix consisting of the coefficients of the highest degree term in 
( /- -1 

each column is nonsingular, i.e., (BG) is nonsingular* But it has 
just been shown that the highest degree of any polynomial in the ith 
column of P(s) is o^. Also from Eq. (3.8), P(s) is given as follov/s 

P(s) = + P 2 ® ■*' ^3®^ 

Thus, from the above representation of P(s), it is clear that the column 
vector consisting of coefficient of S in the ith column of P(s) is 
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nothing but 
represented 
matrix 


the ith column of the component matrix P , ^ which is 


by 


Therefore, P(s) is column proper if the 


iDxra 


(BG) 


1 2 
^^ 0^+1 ^ 02+1 


m 

’'’m+l. 


C4.6) 


is nonsingular. 

Since column properness of P(s) is one of the necessary and 
sufficient conditions for model matching . CSee Sec. 2.5), in summary 
the constraints on the form of P(s) can be stated as follows: 


ti) p 


ai+l 


must be an independent vector for all 


1 < i < m 


(2) P 


i 

c^i+j 


= 0 


for j ^ 2 and for all 1 < i < m 


(4.7) 


The above two constraints in turn limit the realizable set of 
K^’s and require certain sufficiency conditions to be satisfied. These 
conditions are taken into account by Theorem (4,1). The following 
definition is introduced to aid in a precise statement of that theorem. 
Definition 4.1 Component definite polynomial matrix: Let 

(1) the polynomial matrix P(s) be represented by Eq. (3.8), i.e., 

PCs) - Pj^ + + ... (3.8) 


(2) y. be the number of zero columns in P. 
1 1 



68 


(3) r. be a Y, X y. constant matrix formed bv the elements 
1 i 1 

of at the cross section of the v. ordered rows and 
1 1 ■ 

the same ordered columns, which correspond to y^ zero 
columns of 

Then P(s) is called "component definite" if are nonsingular 
for all i >_ 1. 

Notice that to test for component definiteness of P(s) , one need not 

form r. for i = 1 and 2. This is true because each of the a’s is 
1 

greater than or equal to one; therefore, P^ and P^ do not have any 

zero column. Also, since P(s) consists of finite degree polynomials 

only, there exists a finite integer B n, the order of multivariable 

system under discussion, such that P = 0 and thus = P, . Hence, 

p pi 

P(s) is component definite only if P^ is nonsingular. The nonsingu- 
larity of Pj^ is established in Theorem 3.1 and, hence, does not cause 
any new constraint. 

Theorem 4.1 is concerned with the freedom of choice of the 
transfer function elements of the overall transfer matrix T(s). This 
freedom of choice is expressed in terras of the freedom of choice for 
the error coefficients. 

Theorem 4.1 The first a.+l error coefficients of each 

1 

transfer function element t..(s) in ith column of T(s) 

. ji 

can be realized arbitrarily by state variable feedback 
alone if 1) the first error coefficients are chosen 
such that P(s) is forced to be component definite; 2) the 
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last error coefficients , are chosen such that result- 

Oi+l 

ing - are mutually independent vectors for all 1 < m. 

The above theorem gives only sufficient conditions for the freedom 
of choice of the error coefficients of each transfer function. The 
application of the theorem for design of multivariable systems puts 
constraints on only first and last generalized error coefficient for 
each transfer function. The first of the two conditions imposed by 
the theorem is that be chosen such that the resulting PCs) is com- 
ponent definite. The test for component definiteness of P(s) does not 
require lengthy computation because the form of P(s) is completely 
known in advance from the knowledge of the o^'s, and the matrix P^ is 
completely determined from the knowledge of the matrices and 
Thus, one can easily form the x matrices (Def. 4.1) and test 
them for nonsingularity. If some happens to be singular for the 
particular choice of the matrix K^, the must be changed accordingly. 
The second condition of the theorem restricts the last specifiable 
error coefficient to produce a P(s) which is column proper and thus 
realizable. 

The proof is accomplished by first assuming that the conditions 
of the theorem are satisfied, and then showing that all the components 
of P(s) are uniquely determined and that .hey form a P(s) which is 
realizable by the state variable feedback of Eq. (2.2). 
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4,3,1 Proof of Theorem 4,1 

Examine the first equation of the series of Eq. (3.20), 


namely 


KiP^ = 


From Theorem (3.1) it is seen that and must be nonsingular for any 
stable system whose transfer matrix is nonsingular. Thus, knowing 
from the knowledge of R(s) and as specified, P^ is readily given as 



It is assumed that the P^ as obtained above is such as to force P(s) to 
be component definite,. If not, one must change to force P(s) to be 
component definite. 

Now, examine Eq. (3.19) as given below 


K, P . = Q + K P. 

1 r+1 r+1 1 


(3.19) 


where, as denoted also in Eq. (3.19) 


r r+1 


for r = 1 


r-1 

Q = ^ + I ^ .1 . for r > 2 

r r+1 r+l-i i+1 — 


For convenience, let j = r + 1, then Eq. (3.19) can be written as 
follows : 


K, P. = 
1 J 


K.P. 
J 1 


(4.S) 
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Thus in Eq, (4.8), is completely known at the time is being evalu- 
ated, provided is evaluated in strict sequence for j =2, 3, 4, ... . 
This sequential evaluation of the matrices P^ is a significant computa- 
tional advantage. 

Now, examine the form of P(s) . Since the degree of the highest 

degree term in the ith colinan of P(s) is the ith column of P^ is 

identically zero for j > a. +1. Thus, in general, P, might have some 

1 J 

columns identically zero. Let have zero columns, then one need 
only to show that the remaining columns of the matrix are still 

arbitrarily specifiable, because the theorem promises only + 1 arbi- 

trary error coefficients for each element in the ith column of the 
transfer matrix T(s). 

To prove the above, it suffices to show that if i columns 

of are arbitrarily specified, then in Eq. (4.8), the corresponding i 
columns of are uniquely determined from the knowledge of remaining 
m-£ columns of . P^ can then be determined completely from the same 
Eq, (4.8). (Notice that if £ columns of P^ can be specified arbitrarily, 
they could be specified as zero.) The above is proved next. 

Let a^th, a^th, ... columns of P^ be specified arbitrarily 

and let 9^th, O^th, ... 6^ ^th columns of P^ be the remaining unspecified 
columns to be determined. Then by considering the £ specified columns, 
namely a^^th, *** °^£^^ following relationship is 


Kj^[p“^p“2 


p“«]j = 


r ai ao 
[q -Lq ^ 


a (7 T 

q 


K.[p“lp“2 


Ci p 
p 


obtained 
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The above could be written as follows: 


K,W = Z 


where 


W = [p“l p“2 



is a CmxJl) matrix and 


Z = ^ 


0.2 


q""]. 


is a (mxil) matrix. 

ct * 

Since are specified and and are completely knox-m in 

advance, Z is completely known. By expanding the left side of the 
relation K.W = Z, it could be written as follows: 





Z 


where w^ is (lx£) row vector. 

Applying the identity of Eq. (4.5), the left side of the above 
equation could be separated into two parts, one containing the £. un- 
known columns of and the other containing m-Jl arbitrarily specified 
columns of K, , as follows : 
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... k“«], 
J 

^ai 


'^0X 




^02 


w 

otjl ^ 


9m— £ 


(A. 10) 


Now, examine the following IxZ matrix 


w, 


«1 


w 


«2 


w, 




Notice that the matrix W itself consists of Z columns of corresponding 
to zero columns of P , and is the a* the row of W, Thus, by the 

j X i 

definition of component definiteness of P(s) , the matrix under con- 
sideration is nothing but the matrix formed while testing P(s) for a 
component definiteness. Therefore, Eq. (4,10) can be written as follows: 


r. = Z - [k®l k®2 r 


3 3 


w. 


w, 


92 


W0 


m-ilj 


(4.11) 
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Since is an ixZ (here H nonsingular matrix for all j , 

the 9, unspecified columns of on the left side of Eq. (4.11) can be 

determined by multiplying on both sides of Eq. (4.11) by (I'j) 

The matrix is thus completely knoiro, and hence can be 

completely calculated using Eq. (4.8), since everything except is 

known in that equation and is a nonsingular matrix. The matrix 

so obtained is guaranteed to have I columns as specified (identically 

zero in this particular case) , 

Proceeding sequentially for j = 2,3, all the components of 

P(s) can be calculated easily. The P(s) thus obtained is guaranteed to 

have required form (the degree of the highest degree term in the ith 

column of ?(s) is » The only remaining requirement for P(s) to be 

realizable is that it be column proper, i.e., be an independent 

i 

vector for i = 1,2 ... m. However, for some particular choice of 

k^ p^ may not come out to be an independent vector. In that 
i i 

case, one must change k^ to be an independent vector. 

i 

To see that p^ can be made independent by the proper 

i i ^ 

choice of k , - , if suffices to prove that k ,, can be found for any 
a.+l a.+l ^ 

i . X 

arbitrary choice of p But this follows from the fact that Pfs) 

a ,f i 

X 

is component definite. Thus, if while solving for p the columns 

o ^+1 

come out to be dependent, then choose them as desired and solve 

O , "TX 
X 

Eq. (4.8) as before, except that this time Z = Thus, proceeding 

as before, Eq. (4.8) can now be solved for Z = unknown columns of 

J+1 

instead of columns of , The solution is once again guaranteed 
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because F.,- is also nonsingular by the definition of component definite- 
J+1 

ness of P(s). The theorem is thus proved. 

Notice that the constraint that P(s) be component definite 

is easy to handle because this simply requires that P^^ satisfy certain 

properties. However, P^ = and is calculated at the very 

beginning of design. Hence, it is easy to alter. Also, the constraint 

that be chosen so as to force each i to be an independent 

a^+1 *^a^+l ^ 

vector follows from the column properness of P(s). Thus, this con- 
straint amounts to choosing the last error coefficient for each transfer 
function in such a way that- the matrix of Eq. (4.6) as given below 

1 2 m 

Pfli+l Poi+l • ■ • P<J„+1 

is a nonsingular matrix. This requires relatively small changes in the 
error coefficients k 

o.+l. 

1 

In the next section the steps involved in applying the design 
method of this chapter are summarized. 

4.4 Step-by-Step Design Procedure 
The folloxi7ing step-by-step design procedure simplifies the 
presentation of the computational techniques involved in designing the 
multivariable system by using the design method of the chapter. 

Step 1: find the transforma Lj.on matrix Q as follows: 

(a) Find the first n independent columns of the 
controllability matrix. 
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(b) Rearrange these n independent columns to form the 
lexographic matrix and find cf’s from there, 

(c) Take the inverse of this lexographic matrix and form 
the transformation matrix Q as explained in Sec. 2.3. 

Step 2: Apply the transformation to find the new state variable 

representation of the system; that is, find A, B, and C. 

Step 3: Use the results of Step 1 and Step 2 to form the matrices 

S^, S(s), A and B, 

Step 4: Calculate the matrix R(s) and the most general form of the 

matrix P(s) . 

If a decoupled system is desired, follow Steps (5d) through (9d); other- 
wise follow Steps (5) through (12) . 

Step 5d: Test if the system is decouplable and find the pole zero 

excess of each subsystem during this test. 

Step 6d : Find the fixed zeros for each subsystem from the knowledge 

of the matrix R(s). 

Step 7d: Detemine the position of uncontrollable poles, if any, by 

examining R(s) . 

Step 8d: Choose the model transfer matrix for the system which has 

the fixed zeros as found in Step (6d) and pole zero excess 
as found in Step (,5d) , 

Step 9d : Find the pair {F,G> to realize the above decoupled transfer 


matrix. 
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For the general case of model matching, follow the steps listed below: 

Step 5: Expand each t^^(s) of the model transfer matrix in Maclaurin 

series up to a^+1 terms. The specifiable portion of k^'s is 
then completely known. 

Step 6: Write R(s) in its component form. Since R(s) is completely 

known from Step (4) all components of R(s) are known. 

Step 7: Calculate from the knowledge of R^ and K^. 

Step 8: Test PCs) for component definiteness. To do this, form the 

from the knowledge of the form of P(s) and P^. If all 

of them are nonsingular, proceed to Step (9). Otherwise 
change k^ accordingly and go back to Step (7) . 

Step 9: Calculate P^, ... in sequence until all the components of 

P(s) are known. At the end of each sequence, check if the 
p^ are independent vectors. If not, change the 
accordingly and repeat the last sequence. 

Step 10: Find the transfer matrix T(s) = R(s) P ^ (s) , 

Step 11: Find the pair {F,G} to realize the P(s) obtained above. 

Step 12: Using the transformation matrix Q, find the pair {F,G} 

which realizes the transfer matrix of Step (10) , 


The step-by-step design procedure is illustrated with a simple example 
in the next section. 

4.5 Example 4.1 


The example here has been intentionally chosen to be simple in . 
order to meaningfully and concisely illustrate all the steps involved in 
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the design of multivariable systems using generalized error coefficients. 
The example is a slightly changed version of Example 3.1 given in 
Silvinsky’s dissertation [1969] . The problem of decoupling is illus- 
trated first and the problem of complete model matching is tackled next. 

Consider the multivariable system whose block diagram is given 
in Fig, 4.1 and which is described by the following state equations. 


X 


Y 



■ -5 

0 

0 " 


^ 1 

0" 


0 

-1 

0 

X + 

2 

0 


0 

0 

-1 ^ 


^ 0 

1 ^ 


I'l 

0 

01 





X 


L 0 1 1 


u 


(4.11) 


Obviously 

n = order of system = 3 

m - number of inputs = number of outputs - 2 
3 

X e E is a column vector 
2 

Y and U £ E are column vectors 


By inspection, both the input matrix B and output matrix C of Eq. (2.1) 
have full rank, m=2. Also, system is completely controllable and com- 
pletely observable. 

Step 1: The controllability matrix of Eq. (2.7) is obtained as follows: 
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1 0 -5 0 25 0 

2 0 -2 0 2 0 

0 1 0 -2 0 4 J (4.12) 

Clearly, the first n independent columns of the controllability matrix 
are the first three columns of the matrix of Eq. (4,12). 

After rearranging these first three independent columns accord- 
ing to Eq. (2.8), one obtains 

[b^ Ab^ b^] = [b^ (4.13) 

Accordingly, = 2 and = 1, Notice that for this particular example, 

2 2 2 2 

the G^s are unique because the vectors b , Ab , and A b are dependent 
vectors. Hence, = 1, leaving no choice for except that = 2 as 
stated above. 

The inverse of the matrix in Eq . (4.13) is 

0.25 0.625 0 

0,25 0.125 0 

0 0 1.0 J ' (4.14) 

From Eq. (2. 9) , = 2 , = G^ + G^ = 3 and, hence, and are 

defined as the second and third row, respectively, of the matrix in 
Eq. (4.14). 
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Hence, from Eq. (2.10), the transformation matrix is given by 

^1 1 0.125 o' 

= 1.25 -0.125 0 

^2 1 L 0 l.oj (4.15) 

A 

Step 2: Using the tranformation X = QX, the transformed system of Eq. 

(2.11) is found, where the transformed matrices are defined by Eq . 

(2.12) . Thus, the transformed system is given by 




(4.16) 


(4.17) 


As can be checked by inspection, the transformed system matrices 
in Eq. (4.16) do have the special form described in Eq. (2.15), (2.16), 
and (2.17). 

Step 3: The matrices S^, S(s), A and B, are found using Equations (2.18) 

through (2.21), respectively, as given below: 

r s^i 0 1 r o' 
s = = ! 

_ 0 s^2 J Q s J (4.18) 
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S(s) ~ 


1 0 

s 0 

0 1 


(4.19) 


A = 


-5 -6 0 

0 0-2 


(4.20) 


B ^ 


1 0 
0 1 


(4.21) 


Step 4: From Eq. (2,35) R(s) = C S(s). Thus, substituting for C and 


S(s) from Equations 

(4.17) and 

(4.19), 

respectively 


1 

1 0 ' 


' 1 0 ] 

r s-fl 0 ~ 

R(s) = 




1 

= 1 


10 

2 1 _ 


s 0 

L 2(s+5) 1 





- 0 1 . 



(4.22) 


Now, from Eq. (2.36) 


P(s) = (BG)"^ {s'’ - (A + BF) S(s)} 


(A. 23) 


Let 


(BG)"^ = D = 


‘^ll ^12 


^21 ^22 


(4.24) 
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and 


(BG) CA + BF) = -H = - 


^11 ^2 ^13 


^21 ^22 ^23 


(4.25) 


From Equations (4.23), (4.24), and (4.25), one obtains 


P(s) = 


^ '’12® + ‘'ll '‘ 12 ® + ‘'13 


daiS + j 


(4.26) 


Thus, the transfer matrix of the multivariable system in its most 
general form as given by Eq. (2.34) is as follows; 


T(s) = 


s+1 0 


2Cs+5) 1 


-1 


diiS + + ‘'13 


L ‘‘ 21 " ^ ‘' 22 ^ ‘'21 ^22^ + ‘'23 


(4.27) 


4.5.1 The Problem of Decoupling 

In Section 2.6, it was shown that if only the decoupled transfer 
matrix model is desired, then, all the relevant information is obtained 
from the knowledge of R(s) , the system feedback invariant matrix, and 
the necessary and sufficient condition for decoupling of the system 
[Falb and Wolowich, 1967; Gilbert, 1969], 
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Step 5d: If the above mentioned test for decoupling is applied to the 

system in Eq. (4.11) or the transformed system in Eq, (4,16) and (4.17), 
the following information is obtained: 

a) The system is decouplable. 

b) The pole zero excess for both subsystems is 1, 

Step 6d: Now if each row of the matrix R(s) as obtained in Eq. (4.22) 
is written as in Eq, (2.54), the following is obtained: 

R^(s) = [s+1 0] = (s+l)[l 0] 

R 2 (s) = [2(s+5) 1] = l[2(s+5) 1] (4.28) 

Thus (s+1) and 1 are the highest degree polynomials common to all entries 
of the first and second rows of R(s) respectively. Since the zeros of 
the highest degree common polynomial of ith row of R(s) are the zeros 
of ith subsystem the following information is obtained: 

a) The zeros of the first subsystem are the zeros of s+1. 

b) The second subsystem has no zeros . 

Once again, since the order of ith subsystem is the sum of the number of 
zeros and pole zero excess, the following conclusion can be drawn easily. 

a) The order of the first subsystem is 1+1 = 2, 

b) The order of the second subsystem is 0+1 = 1. 

Step 7d: ■ Since the sum of the orders of the two subsystems = 2+1 = 3 
is equal to the order of the overall system of Eq. (4.11), there are no 
uncontrollable poles. This fact could also be verified by checking the 
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matrix R’ (s) of Eq. (2.62). The uncontrollable poles are the zeros of 
(r^ (s) I . For this example 

|r’(s)| = 

Since the polynomial |r’Cs)| has no zeros, there are no uncontrollable 
poles as stated before. 

Step 8d: In summary, if decoupling is desired then choose the model 

transfer matrix such that the first subsystem has a zero at s = -1 and 
two arbitrary poles and the second subsystem has no zero but one arbi- 
trary pole. 

Step 9d: Hence, given a transfer matrix model which meets the above 

specif ications , the pair {F,G}, to realize the model response from the 
plant, could easily be found by using any one of the following methods: 

1. Original algorithm of Gilbert [1969] for decoupling 

2. Wolowich^s [1972] algorithm for exact model matching 

3. Design method desicribed in this study. 

All of the above algorithms give the same F arid G, Method 2 is compu- 
tationally more difficult because, as mentioned in Sec. 2.5, it requires 
that the model transfer matrix T(s)' be put into the form 

T (s) = R (s) P ”\s), 

m mm 

where and are relatively right prime polynomial matrices as 


2 (s+5) 


= 1 


given in Eq . (2.38). 
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The algorithm for finding the feedback and gain matrices F and 
G is illustrated for the general case of model matching which is con- 
sidered next. 


4.5,2 The General Case of Model Matching 

The general case of model matching is illustrated by trying to 
approximately match the plant transfer matrix to a given model transfer 
matrix using state variable feedback. Assume that the following model 
meets the design requirements for the given plant of Eq. (4.11). 


T(s) 


6(s+l) 
(s+2) (s+3) 

10 

^ (s+2) (s+5) 


0.6(s+l) 
(s+2) (s+3) 

4 

(s+4) 


(4,29) 


If one were to test for exact model matching, it would be found that the 
necessary and sufficient conditions for exact model matching as given in 
Theorem 2,2 IWolowich, 1972] are not satisfied, and, hence, the plant of 
Eq. (4.11) cannot be matched to the above model. Also, application of 
the test does not give any hint as to what changes should be made in the 
model to force realizability. Thus, the designer is left to his luck 
and experience to try one model after another until he finds one that can 
be matched and that meets his requirements. 

The design method described in the previous settion is now 
applied to investigate if there is any otl-er model which matches the 
plant of Eq. (4.11) and at the same time, approximates the model of 
Eq. (4.29). 


\ 





G 






Not 
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Step 7: Substituting for and from Eq. (4.30) and (4,31), one 

obtains by using Eq. (3.20) 




Hence, 


^1 = 


' 1.0 

0.1 ^ 

“1 

“ 1 

0 “ 

1.0 

1.0 _ 


_10 

1 _ 


0 

10.0 


-.111 

1,111 


(4.32) 


Step 8: To check for the component definiteness of the polynomial 

matrix P(s) j expand P(s) as obtained in Eq . (4.26) in its components 
as described in Eq. (3.8). Clearly, 

P(s) = + P^s + 


Since the are greater than or equal to 1, P^ and P^ do not have 

any zero columns, and hence, one need not form and Since the 

last non-zero component of F(s) is P^, one needs to form and check 
only. From Eq. (4.26) 



0 


(4.33) 
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and, hence, only the second column of is zero. Thus, according to 
definition 4.1, is a 1x1 submatrix formed by the elements at the 
cross section of the second row and the second column of Thus, 

Fj = (P22)i =1.111 


and, hence, is nonsingular and thus P(s) is component definite 
matrix. 

Step 9: Now find using Eq. (3.20) 

V 2 = '^2 + ‘^2 h . 


or 


P2 = {R2 + 


K. 




After substituting for K^, P^, and R 2 , the result is 


P = 
2 


1.0 0.1 


‘0 
^ 1— 1 


-0.167 -.0167 


' 0 -.111 


i 


+ 1 




1.0 1.0_ 

\ 

2 0_ 


0.70 0.250 


10.0 1.111 


Hence, from above 



^1.0 

0.1 

-1 

0.833 

a ' 

•”2 “ 

1.0 

1.0 


2.50 

0.2 


0.426 -.022 


4.07 


0.222 


(4.34) 
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Now, from Eq. (4.6), P(s) is column proper if 


1 2 
^0^+1 ^0^+1 


contains all independent vectors, i.e., the columns of the following 
matrix 


1 2 
P3 P2 


(4.35) 


2 

must be Independent vectors. Since is known, one must check that it 

be a non-zero vector. From Eq. (4.34) it is seen to be a non*~zero 
vector; hence, proceed to find . 

Solving Eq. (4.8) for r=2, one obtains 


P3 = Q2 + K3P^ 


(4.36) 


and 


'^2 = ^"3 + ’'2 ^2 


Since R^, and P 2 are all known at this point, Q 2 is obtained by direct 
substitution as fellows: 


Q2 = lOJ + 


i.e. , 


Q 2 - 


-.1667 -.01671 

0.70 0.250, 

' -.139 0 

^ 1.32 0.04, 


0.426 -.022 

4.07 0.222 
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Next, from Eq. (4.30), since the second column of is unspeci- 
fiable, this must be determined by considering the second column of both 
sides of Eq. C4-36). 



2 

’2 


+ K 


2 

3P1 


But, from Eq. (4.31), p = 0 and hence 


2 

'3 '*1 


K., = - q 


or 




k.- 


'"(Pi2h’ 

^^22^ 


= - q< 


The above can be partitioned into known and unknown parts as shown in 
Eq. (4.10). 

S ^^12^1**' S ^^22^1 " ^ ^2 

Hence 
and thus 

kj = “|92 ^^12^lJ [^^22^1] 
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Substituting for all quantities on the right hand side of the above ex- 
pression, one obtains 


k 


2 

3 



~ 0 ^ 


0.306“ 

1 



+ 


[-0.111] 

1 

_0.4_ 


-.390 

J 


[ 1 . 111 ] 


-.0340 

0.0833 


( 1 . 11 ) 


0.0306 

-.0750 


The remaining first column of (the only unknown part) is now 
calculated easily from Eq . C4.36) by either solving for the first column 
of or by solving for the complete P^, even though the latter involves 
unnecessary more computation. So consider the first column of both 
sides of Eq. (4.36). 


V 1 1 1 

P3 = q2 + K3 


-0.139“ 

1.32 


0.306 
, -.39 

0.0306“ 

-.0750_ 

“ o' 
10 

-0.139' 

1.32 

+ 

0.306“ 

-.750_ 

= 

“0.167 

0.567_ 




ZJ 
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Hence 


1 

^3 



0,1 67~ 
0.667 


1.0 

O.ll 

-1 

0,167“ 

1.0 

i.qj 


0.567: 


0.123^ 

0.444 


Now, from Eq. (4.35), must be independent of p^. Since p^ is 

2 1 

the second column of in Eq, (4,34), by inspection and are 
mutually independent column vectors. 

Except for the determination- of the necessary F and G, the de~ 
sign is thus complete, since PCs) is completely known and is guaranteed 
to satisfy the generalized error coefficients K^, K^, and as speci- 
fied in Eq . (4.30). Next, check the designed transfer matrix so 
obtained . 

Step 10; The polynomial matrix P(s) is completely specified as follows; 


PCs) = ^ 3 ®^ 


,123s + ,426s' 


-C.0222S + .111) 


.444s^ + 4,07s + 10 .222s + 1,111 


(4.37) 


Hence 


P ^(s) = P*(s)/A 
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where 


P*(s) = adj P(s) 


\?.22(s+5) ,0222(s+5)‘ 

_(.444s^ + 4.07s + 10 ,123s^ + .426s 


(4.38) 


A = |P(s)| = (.123s^ + .426s) (. 222) (s+5) + (.444s^ + 4.07s + 10) ( . 022) (s+5) 
= .022 (s+5) {1.23s + 4.26s + .444s + 4.07s + 10} 

= .022 (s+5) U.674 s^ + 8.33s + 10} 

= .0371 (s+5) {s^ + 5s + 6} 

= .0371 (s+5)(s+2)(s+3) 

From Eq. (2.34) 

T(s) = R(s) P“^(s) 

= R(s) P*(s)/A 

Substituting for R(s) and P*(s) from Equations (4.22) and (4.38), 
respectively 

-.222(s+5) .0222(s+5) 

-(.444s^ + 4.07s + 10) .123s^ + .426s_ 

“'.222(s+l) (s+5) .0222(s+l)(s+5) 1 

.444 (s+5) ^ - ,444s^ - 4.07s “ 10 .0444 (s+5) V.123s^+.426sJ 


R(s)P*(s) = 


(s+1) . 0 

2 (s+5) 1 
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“.222 (s+l)(s+5) .0222 (s+1) (s+5) 

.37s + 1.111 .167s^ + .870s + 1.111 


Therefore 


T(s) = R(s) P*(s)/A 


'.222 

(s+1) (s+5) 

.0222 (s+1) (s+5) 

.0371 

(s+5) (s+2) (s+3) 

.0371 (s+5)(s+2)(s+3) 


.37 (s+3) 

.167 (s+2. 22) (s+2. 98) 

.0371 

(s+5) (s+2) (s+3) 

.0371 (s+5) (s+2) (s+3) 


6 (s+1) 
(s+2)(s+3) 


0.6 (s-fl) 
(s+2) (s+3) 


- 10 

[(s+2)Cs+5) 


4.5(s + 2.22) 
(s+5) (s+2) 


(4.39) 


Some comments are in order before the pair {F,G} are found to 
realize the transfer matrix of Eq. (4.39). The transfer matrix actually 
desired is given in Eq. (4.29). A comparison between the transfer 
matrices of Equations (4.29) and (4.39) shows that except for the trans- 
fer function entry ^22^^^ other entries are exactly realized. Also, 

for the transfer function t 22 (.s)» the difference between desired and 
obtained is not intolerable, as can be seen by comparing the time re- 
sponse of the two as shown in Fig. (4.2). 

However, one may not always be as lucky as in this example. As 
mentioned in the beginning, this is a relatively simple example to 
demonstrate the design method of multivariable systems. It is quite 
possible and probable that the response obtained would not be as close 
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to the desired response as indicated in Fig, 4,2 for this case. For 
this example one could guess to some extent which of the transfer func- 
tion entry would differ most in its response. Since the second column 
is not specifiable, one would guess some change in the second 
column of Tfs), Also, since t^ 2 (®) just a constant multiplication 
of tlfCs), one would not expect much change in Thus, the only 

entry left is ^ 22 ^^^ which might have much different transient response 
than desired. As can be seen by examining 

(4.39), they both have closed loop poles at the same position. (Thus, 
it seems (^ 21^3 more effect on ^ 22 ^^^ than any other single error 
coefficients,) The above comments are included to give the reader some 
insight to the problem solution. However, extreme caution should be 
taken before any of the above concepts are generalized. 

Once again, it should be noted that even if the desired transfer 
function and the one obtained might differ significantly in their tran- 
sient behavior, they still asymptotically coincide in their steady state 
behavior if. the system design is stable. 

Step 11: To find the pair {F,G} which realizes the transfer matrix as 

obtained in Eq. (4.39), find the pair {F,G} to realize P(s) as obtained 
in Eq. (4.37). . 

Comparing Eq. (4.26) and (4.37), the D matrix of Eq. (4.24) is 
given as follows: 


pll 

^12 


\l23 

-. 0222 ' 

h^ 2 i 

c. 

^ 22 ^ 


.444 

.222 
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Hence, from Eq. (4.24) 


(BG) = D 


-1 


'.123 

-.022~ 

-1 

6.0 

.60 

.444 

.222_ 


-12.0 

3.3 


(4,40) 


Hence, G=B^D^=D^, because from Eq, (4.21) B = Identity matrix 
Thus 


G = 


6.0 

.60 

[-12.0 

3.3 


(4.41) 


Once again comparing Equations (4.26) and (4.37), the H matrix of 
Eq. (4.25) is given by 


H = 



^12 

h3 


“ 0 

.426 

-.11 


^22 

^23 


10 

4.07 

1.11 


(4.42) 


But, from Eq. (4.25) 


Hence 


This implies 


(BG) ^ (A + BF) = -H 


A + BF = -(BG)H 


F = { (BG)H + A} 


= --GH -BA 
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After substituting for G, H, B and A from Equations (4.41), (4.42), 
(4.21) and (4.20), respectively, one gets 


F = 


^ 6.0 


6 ' 

■ 0 

.426 

-.11 

- 12.0 

3 . 3 J 

10 

L. 

4 , 

07 

1.11 

1 0 ~ 


' -5 

- 

-6 

o “ 


0 1 ^ 


0 


0 

-1 


6.0 5.0 

0 l 


r-5 

-6 


L33.0 8.37 5 


0 0 -2 


- 1.0 1.0 o ' 

- 33.0 - 8.37 - 3.0 


Step 12: Using Eq. (2.14) 

F = FQ 


(4.43) 


(2.14) 


After substituting for F and Q from Equations (4.43) and (4.15) 


respectively, one obtains 
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'-1.0 1.0 o' 


-.25 .125 o' 



1.25 -.125 0 

-33.0 -8.37 -30 




0 0 1 


ri.50 -.25 0 


[-2.2 -3.08 -3.0 


(A. 44) 


Thus, Equations (4.44) and (4.41) give the required feedback pair {F,G} 
which realizes the transfer matrix of Eq. (4.39). 

This satisfactorily completes the design. The state variable 
feedback compensated plant is shown in Fig. (4.3). In the next section 
the results of this chapter are summarized, and advantages and disadvan- 
tages are pointed out. 

4.6 Conclusion 

A completely new design method for the design of multivariable 
systems is presented. The main advantage of the above design method is 
that it realizes a given (pre-specified) model transfer matrix approxi- 
mate, where necessary and sufficient conditions for exact model matching 
fail to produce any solution. The design method is based upon approxi- 
mating each element of the transfer matrix by a finite element series. 

The coefficients of the elements in the series are well known to have 
direct relationship with the generalized error coefficients. As compared 
to exact model matching, this design method tries to realize generalized 
error coefficients of the model. It is shown that certain sufficiency 
conditions (Theorem 4.1) must be met which in turn restrict the allowable 
















range of error coefficients. This in itself is not a serious design 
restriction for all practical purposes. 

The main disadvantage of the design method is that the design 
does not take into account the transient behavior of the system. This 
problem could theoretically be serious enough to cause instability. 
This disadvantage is attributed to the fact that an n terms series 
representation of an nth order transfer function element does not 
uniquely determine the transfer function. That is to say, there are 
many transfer functions, some of them having undesirable transient 
behavior or, in the worst case, representing unstable systems, which 
have the same n generalized error coefficients. 

Thus, it is possible that the application of the above design 
procedure may result in an undesirable transient response. Of course 
this fact may be checked easily by comparing the desired and actual 
step function responses. In the next chapter a second design method is 
given which ensures stability on the risk of deteriorating the steady 
state behavior of some of the elements in transfer matrix. Thus, even 
though an overall satisfactory solution is not guaranteed in advance, 
the design method of this chapter gives valuable insight into the 
solution to the problem of model matching. 



CHAPTER 5 


THE POLE FIXING METHOD 

5.1 Introduction and Organization of the Chapter 
The pole fixing method for the design of multivariable systems 
is the subject matter of this chapter. The pole fixing method has a 
distinct advantage over the error coefficient design method of the last 
chapter. The shortcomings of the error coefficient method are summa- 
rized and the salient features of the pole fixing methods are discussed 
in the next section. The constraints of the pole fixing method are 
described via Theorem 5M in Section 5.3. The. proof of the theorem is 
constructive and leads to a design procedure. The step-by-step design 
procedure is described in Section 5.4. In Section 5.5 the pole fixing 
design procedure is illustrated by reworking the example of the last' 
chapter. Some special cases of the pole fixing method are discussed in 
the next section. Finally, the findings of the chapter are summarized 
in Section 5.7. 


5.2 Need for Pole Fixing 

The error coefficient design method of the last chapter is based 
completely upon the generalized error coefficient representation of a 
transfer function. The generalized error coefficients describe the sys- 
tem input behavior only after sufficient time has elapsed for the system 
transients to decay to insignificant amplitudes. Thus, even though the 
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error coefficient design method of the last chapter is simple, it does 
not ensure satisfactory transient behavior of the systems. In the worst 
case the design may even lead to an unstable system. 

The pole fixing method to be introduced in this chapter ensures 
the system stability by prefixing all the system poles at specified 
locations in the left hand S plane. It is well known [Anderson and Luen- 
berger, 1967; Wonham, 1967; Davison, 1968; and Sridhar and Lindorf f , 

1972] that if the multivariable system of Eq. (2.1) is completely con- 
trollable, then all of the n poles of the system can be fixed arbitrar- 
ily by using l.s.v.f. alone; The pole fixing method, in addition to 
assigning all the poles of the system, fixes zeros and steady state errors 
for some of the transfer functions in the transfer matrix. However, the 
fixing of the zeros is done indirectly via error coefficients, and re- 
quires that , the first component of the system invariant matrix R(s) , 
satisfy certain suffficiency conditions. Most of the time these condi- 
tions can be met by simple changes in the names of inputs and outputs. 

Compared to the decoupling technique, the pole fixing method has 
three distinct advantages: 

1, considerably more freedom for the off-diagonal terms 
is achieved, 

2, controllable . systems which cannot be decoupled can be 
designed for undecoupled response, and 

3, there are no uncontrollable poles. 

Compared to the error coefficient method of the last chapter, it has the 
advantage that the transient behavior can be controlled as well. The 
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only disadvantage is that no error coefficients can be specified for one 
side of the diagonal terms In the transfer function matrix which amounts 
to not being able to control the zeros of that side of the transfer 
function matrix, directly or indirectly. However, if the open loop 
system has only one way coupling, i.e., R(s) is a triangular matrix, or 
can be made a triangular matrix by changing the names of inputs and out- 
puts, then one side of the diagonal terms in the transfer matrix could 
be made identically zero. This is discussed in Section 5.6 under the 
special case of triangular decoupling - 

Thus, despite the restrictions that the zeros for only one side 
of the diagonal terms in the transfer function matrix can be specified, 
the pole fixing design method is useful for two reasons. First, it gives 
a lead into how the transient behavior specifications can be accoiranodated 
in the design criteria, and, secondly, information is obtained as to 
what is realizable from the system. 

In the next section the constraints of the pole fixing method are 
mentioned in Theorem 5.1. The proof of the theorem leads to a design 
procedure. 

5.3 The Pole Fixing Method 

In this section the design constraints are introduced through 
Theorem 5.1. The design procedure follows from the proof of the theorem. 
First, a definition is introducer! to help understand the implications of 
the theorem that follows: 
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Definition 5.1: Lower (upper) Definite Matrix: An mxm 

constant matrix P is called Lower (upper) definite if the 
YXY submatrix formed from the lower right most (upper 
left most) terms of P is nonsingular for all 1 ^ Y ^ 

A matrix which is both lower definite and upper definite is simply 
called definite nonsingular matrix . Conversely, if a matrix is 
definite nonsingular, it is both lower definite and upper definite. 
Notice that any triangular matrix (one side of diagonal terms iden- 
tically zero) if nonsingular is a definite nonsingular matrix, a 
requirement which is satisfied if, and only if, none of the diagonal 
terms are zero. 

Theorem 5.1; If the multivariable system is controllable, 

then all the n poles of the overall system can be placed 

arbitrarily via l.s.v.f. along by forcing P(s) in Eq . (2.36) 

to be triangular. Moreover, if R, = R(s) is lower 

i s=o 

(upper) definite, where R(s) is the system invariant matrix 
of Eq. (2.35), then the following freedom in the choice of 
error coefficient is maintained by choosing P(s) lower (upper) 
triangular . 

(1) Elements of can be chosen arbitrarily, for i ^ j (i < j), 

provided the elements (k.,)- of K_ are so chosen as to force 

xj 1 1 

to be a lower (upper) definite matrix. 

(2) addition error coefficients can be arbitrarily assigned 
for those elements t..(s) of the overall transfer matrix T(s) 
for which i>j(i<j). 
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Throughout the statement of the above theorem, the expression in paren- 
thesis corresponds to being an upper definite matrix. Notice that if 
is definite nonsingular, (i.e., is both lower and upper nonsingular) 
then P(s) could be chosen either upper triangular or lower triangular. 
This results in additional design freedom, as illustrated in the example 
of Section 5.5. 

Also notice that since P(s) is forced to be triangular and since 

the ith diagonal element in P(s) is chosen to be a polynomial of degree 

Cessential for choosing all n poles), P(s) Is automatically forced to 

be column proper. Thus, P(s) as specified above is realizable by linear 

state variable feedback alone. 

Proof of Theorem 5.1: The first part of the theorem is proved 

by noticing that if the system is controllable, then the 
Lunenberger transformation of Section 2,3 can be found and the 
overall transfer matrix of the system can be written as 
T(s) = R(s) P ^(s) as given by Eq. (2,34). Thus, all n poles 
can be placed arbitrarily simply by forcing P(s) triangular 
with ith diagonal term a polynomial of degree where the 
<7^ are specified by Eq, (2.8). But such a P(s) is realizable 
by l.s.v.f. alone because it is column proper and ith column 
has the ith diagonal term as the highest degree polynomial in 
it. Thus, the l.s.v.f, pair (F,G) can be found by simply 
equating the desired P(s) with the one in Eq, (2.36). 

The second part of the theorem is proved by showing that 
the off diagonal nonzero terms in P(s) can be used to specify 



108 


the error coefficients of the corresponding terms in the 
transfer matrix. Conversely, it is shown that the off 
diagonal nonzero terms in P(s) are calculated for any 
arbitrary choice of the error coefficients as specified in the 
theorem. This is accomplished in two parts. First, the matrix 
the first component of P(s), is determined. Notice that 
to specify all n poles, the diagonal terms in the triangular 
matrix P(s) are specified to within a constant only. This 
gives extra freedom, i.e., the first error coefficient can be 
specified even for the diagonal terms in T(s) . Thus, P^ is 
determined first. Next, the remaining components of P(s) are 
determined one by one. Notice that only terms on one side of 
the diagonal need be determined because the diagonal terms are 
now completely specified and the terms on the other side of 
diagonal terms are specified to be identically zero. 

The theorem is proved for an upper definite matrix. The 

proof for lower definite can be developed on the same lines. For R^ 
upper definite, P(s) is forced to be upper triangular and, hence, , 

P^, ... are all upper triangular matrices. 

To determine P^ , consider Eq, (3.19) for r = o, which is as 

follows : 

K^P^= (3.19) 

where P^ and are ail mxm constant matrices. 
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Now, consider the column of both sides of the above equa- 
tion, Since is an upper triangular matrix, the ith column is 
simply 


ir 


k”j 


li 


= [r^]. 



The above could be simplified as follows, where zero terms have been 
dropped. - 


[k^ ... k^] 


Pll 


P ■ ♦ 

11 


1 - 


= [r,] 


Next, consider the top i rows and the. remaining m-i rows of the above 
equation separately as follows: 


11 

... k, . 

li 


Pii 


’^li 

• 



• 


* 

^il 

... k. , 

ii_ 

1 

J’ii. 

1 

r . . 

ii_ 


( 5 . 1 ) 
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and 


^i+1 1 

\+l i 


Pli 


^i+li 

1 

... k 

m X 


Pii 


r 

m 1 

- 

- 

1 

J 

1 

_ 


(5.2) 


Now, Equations (5.1) and (5.2) must be solved for i = 1, 2 ... m. 
Substitution of i = 1 in Equations (5.1) and (5.2), respectively gives 


(ki^)l (PiPi = 


and 


^21 

Ir 


^21 

ml 

1 

/ml 


Since R (s) and, hence, are already known, (p-j^-j^)j^ 


^12 


'Im 


are 


are uniquely determined for any non 2 ero choice of 

Equation (5.1) can next be solved -for i == 2 to give 



^12 

■ 

pi 2 ' 


’^12 

/ 2 I 

^22_ 

1 

|/22_ 

1 

r^2 


from which 


'12 


22 


L 


is uniquely determined for any choice of 



Ill 


^^22^1 keeps 


•^ll 

CM 

‘‘21 

^"22 


'-■1 


a nonsingular matrix. Notice that since 


. Hence, by 


is upper definite matrix 

— - 
^11 

is independent of | 

^12 


^12 

i 

1 i 

^^22 


considering the first two rows and top two columns of Eq. (3.19), it is 


easily seen that 




^12 

1 

is independent of 

r “1 
»ii 

^22J 

1 

0 


which in turn implies 


If m > 2, one would need to solve Eq . (5.2) for the remaining 
unspecif iable elements in the second column of k^. This can be done by 
splitting the left side of Eq, (5.2) into two parts as follows 


^i+1 1 

• » • k , , _ 

1+1 

i-1 


Pii 


'“i+i i‘ 


^i+l,i 

* 




• 

+ 

: 

• 

^Piih = 

* 

k _ 
ml 

... k . 

m 1 - 

1 

1 

^i~l i 

1 


1 

\ i 


and the above can be rewritten as follows 


^i+1 , i 


_ 1 


^i+1 i 


^i+1 1 

— 

•••. ^±+1 i_i 


Pi 1 ’ 


% 

k 

m i _ 

1 

(pii>i 

L 

r 

mi 

1 

k 

m, 1 

k . 1 
m 1-1 

1 


> 

-1 1 


(5.3) 


Also, Eq. (5.1) can be re^^ritten as follows 
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11 


il 



(5*4) 


Thus solving Eq. (5.3) for i = 2, gives the remaining elements of the 
second column of . 

The above result can be used to solve Eq. (5.4) for i = 3, and 
the results could then be used to calculate the remaining elements of 
the third column of if m > 3. This sequence is repeated until Eq. 
(5.4) is solved for i = m which gives the desired upper triangular 
for the desired upper triangular portion of (Including the diagonal 
terms). 

Once again, it can be seen from Eq. (5.4) that a unique solution 
of Eq. (5.4) requires that elements in be chosen such that is an 
upper definite matrix. Similarly, a unique solution of Eq . (5.3) 
requires that (p^^)^ ^ which is obtained only if is upper definite. 
This proves part (1) of Theorem (5.1). 

To prove the second part of the theorem, that is, that the remain- 
ing components of P(s) are uniquely determined by an arbitrary choice 
of corresponding components of T(s) , consider Eq . (3.19) for r > 1, 
which is given below as 


1 r+1 ^r r+1 1 


(3.19) 
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Notice that if the above equation is solved in sequence for 
r=l, 2, ...» then is always known. Also, and are completely 
known from above, and the diagonal terms in completely known 

because of the fact that the diagonal terms in P(s) are now completely 
specified. Finally, is upper triangular because P(s) is assumed 

to be upper triangular. Thus, it suffices to prove that the terms above 
the diagonal term in ^nd terms below and Including diagonal terms 

in are uniquely determined for arbitrarily specified terms above 

the diagonal in 1o see this, let j = r+1 and rewrite the above 

equation as follows: 

K,Pj=Q^ + Yl (5.5) 

Consider the ith column of both sides of the above equation which is 
given as follows: 


[k^ ... k“] 

^li 

= + [k^ ' • • k^] 

r J 3 

Pii 


• 

« 

Pii 


Pii 


0 

■ 

0 


P _ 

J 

0 


(5.6) 


Using identity (4.1), the above equation could be written as follows: 
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[kj . . . kj”^] 

i ' 

+ IkJ] (PjPj + [kj+^ 

... k^] 

1 

^ 0 " 


h-1 i 


] 

0 


= Iqp + 



“r‘> 

Pi i ' 

+ Ik^h (p.Pi + [kj+^ 

...k“] 

o' 





0 


Pi^l i 

1 




C5,7) 


which can be simplified. to give 


IK 




Pi-i 1 


+ [k^] (P,,) . 


[qh+ tkl^ 

**" J 


J 


li 


Pi--1 i 


+[kJj(Pii)i 


(5.8) 


Now, break the above Eq, (5.8) into two parts. 

1. Recombine the last two terms of the right side of the equation 
and consider the top (i-1) rows of the both sides of the 
equation thus obtained. The following results 
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2. Again, recombine the first two terms of the left side of 
Eq. (5.8) and consider the last m-i+1 rows of the result. 
The following equation is obtained. 



(5.10) 
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If [pji is not known completely, then the only unknown part is 
the terms above the diagonal term in that column which could be uniquely 
determined by using Eq. (5.9) for corresponding arbitrarily specified 
terms in This is so because all other terms are known in Eq. (5.9), 

and since being an upper definite matrix, inverse of its upper sub- 
matrix exists. 

Finally, the remaining portion of the ith column 
uniquely determined from Eq. (5.10) because 0 snd all other 

quantities are known. 


Thus, the ith column of and are uniquely determined for 

arbitrary above the diagonal terms in • Incrementing i and repeating 

the procedure until i = m, completely specifies with as stated 
in the theorem. Next, increment r, and solve for the unknown part 


of the next higher component of P(s) until all the components are deter- 
mined. The P(s) so obtained remains column proper with the degree of the 
highest degree polynomial in ith column being a^. Hence, the l.s.v.f . 
pair {F,G} can be found to realize the above P(s), which in turn forces 
all the n poles of the system as desired and achieves certain other 
desired properties as specified in terms of error coefficients. 

The theorem is thus proved for the case of a upper definite 
matrix. For the case of a lower definite matrix one can proceed 
similarly by forcing P(s) to be a lower triangular matrix. The equations 
involved in solving for the components of P(s) are given below without 
proof. 


To calculate and , let 2, = m-i+I and solve the following 


tl XOnS S jlOjT jl “ X*j /L « • * ill* 
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(5.12) 


The use of Equations (5.11) and (5.12) in solving for and unknown 

part of k^ is further demonstrated via example of Section 5.5. 

To solve for the unknown parts of P. and K. , let Jl = m-i+l and 

1 1 

proceed in sequence for i = 1,2 ... m. 

£ 

If at any point [p^] is completely specified, then determine 

a 

[k^] from the following 




(5.13) 
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Otherwise, calculate unspecified parts of p^ for corresponding 

a 

arbitrary elements in as follows; 



k ■ 1 

• * ' M m 

-Ij 

1 



+ 

^£+1 m 

« 



* * * 

m m 

1 


V 
! 

r 

ic A » « « Ic 
m £ mm 

j 

^ml 

L J 


't 


£+1 £ 


• 


£ 

1 


(5.14) 


and , then , 
equation. 


finally find unspecif iable parts of k. using the 


following 


1 

1 



••• ^m 


P££ 


^1£ 


^1£+1’ * '^Im 


^£+1 



• 



• 

• 


• 


• 


• 


• 



^££ 

j 


^££ 

£m 


Pm£ 


^££ 


_^££ ^£m 

j 

£ 



1- ^ 


L. 


1 

U. 

j 


r 






(5.15)' 

In all of the above equations H = m-i+1, and one must proceed sequen- 
tially for i = 1, 2 ... m as chov/n in the example of Section 5.5. 

The proof of the theorem is thus complete. In the next section 
a step-by-step design procedure is outlined which is used in a subse- 


quent example. 
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5.4 Summary of the Design Procedure 
The statement and proof of Theorem 5.1 has resulted in^design 
constraints and design procedure. The constraints are summarized below 
and are followed by step-by-step design procedure for quick reference. 

The main advantage of the design method described in this section 
is that all the n poles of overall system are specified (fixed) in 
advance by forcing P(s) a triangular matrix and by specifying all diag- 
onal terms to within a constant. In addition, if R(s) meets certain 
prespecified conditions, then additional error coefficients can be 
specified for those terms in TCs) for which corresponding terms in P(s) 
are free. Thus, if P(s) is forced to be a lower (upper) triangular 
matrix, then one has the following information: 

1. The terms above (below) the diagonal terms in P(s) are forced 
to be zero and, hence, for corresponding above (below) the 
diagonal terms in T(s) no error coefficients can be specified. 

2. The diagonal terms are fixed to within a constant, and thus 
these constants can be used to specify the first error 
coefficients (position error coefficients) of the diagonal 
terms in T(s) . 

3. The terms below (above) the diagonal terms in P(s) are free 
polynomials of degree _< for the ith column in P(s), and 
hence can be used to specify o^ more error coefficients for 
the terms in ith column of T(s) . 


Thus, for each component of ?(s) except for P^ only the terms below 
(above) the diagonal are unknown and all other terms are known. For 
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all terms below (above) and including the diagonal terms are unknown 
and the rest are identically zero. For all terms below Cabove) and 
including the diagonal terms can be specified as desired provided they 
are specified in such a way as to force lower (upper) definite. For 
, j ^ 2 only terms corresponding to the free elements in can be 
specified and the rest must be determined. - 

The step-by-step procedure below utilizes the above information 
and summarizes the steps involved in the design. 

Steps 1-5; These steps are exactly the same as those given for 
the design method of the last chapter. These consist of deriving the 
feedback invariants of the system. 

Step 6: Break R(s) in its components and examine . If is 

lower Cupper) definite, choose P(s) lower (upper) triangular. If is 

definite nonsingular P(s) could be chosen either lower or upper triangu- 
lar, whichever form is more useful. Go to Step 8. 

Step 7: If is neither upper triangular nor lower triangular, 

then a change in the names of the inputs or outputs might do the job in 
most cases. If so, repeat Steps 1 through 6. 

Step 8: Determine and the unspecif iable portion of by 

using Equations (5.3) and (5. A) as described in the last section, 

[Use Equations (5,11) and (5.12) for lower triangular P(s).] 

Step 9; Determine the diagonal terms of P(s) from the knowledge 
of the diagonal terms in P^. The diagonal terms in P(s) are specified 
to within a constant and these constants can be determined by comparing 
them x^ith the- diagonal terms of P^, 



121 


Step 10: Find the complete and by using Equations (5,8), 

(5,9) and (5.10) and solving for i = 1, 2 m in sequence. [Use 

Equations (5.13), (5.14), and (5.15) for P(s) lower triangular.] Solve 
for j =1, 2 ... until all components of are completely known. 

Step 11: Find the closed loop transfer matrix T(s) = R(s)P ^(s) 
as both R(s) and P(s) are now completely known. 

Step 12: Find the l.s.v.f. pair {F,G} to realize this T(s) using 

Eq. (2,36), 

This completes the design. In the next section the example of the last 
chapter is reworked using the new design method of this chapter. 


5.5 Example 5.1 

Consider the same system as in the example of the last chapter. 
Since Steps 1 through 5 in the pole fixing method are exactly the same . 
as in design method of the last chapter, the results of these steps 
are summarized below. For details, see Section 4.5. 


Step 1: 


= 2, = I 


Step 2: The transformation matrix Q is 


Q = 


- .25 0.125 0 

1.25 -.125 0 

0 0 1.0 
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Step 3: 



0 

-1 

1 

' 1 

o' 


s 

0 




S*^ = 

0 

s 

, S(s) = 

s 

0 


- 

J 


0 

1— 

1 _ 


Step 4: 


R(s) = C S(s) 


’ s+1 0 ~ 

_2(s+5) 1 


°11 ® ^11 ^ 12 *^ ^3 


P(s) = 


d2is + + h^3_ 


^22® ^23 J 


Step 5: Break the elements of the desired transfer matrix in their 
Maclaurin series. This is done in Section 4.5, and is not repeated here 


Step 6: 


- R(s) 


s=o 


s+1 


2(s+5) 1 


1 0 

1C 1 


s~o 


Similarly, 


^2.= 


1 0 "^ 
2 0 


, R3 = = ... = 0 


(5.16) 
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By inspection, R is definite nonsingular matrix and, hence, is 
both lower and upper nonsingular (Def, 5.1). P(s) could thus be forced 
either lower or upper definite. As mentioned earlier, this is an added 
advantage. Since = 2 and =1, forcing PCs) to be lower triangular 
results in greater design freedom because for t^^Cs), + 1 = 3. Hence, 

three error coefficients can be chosen. By -forcing P(s) upper triangular, 
only + 1 = 2 error coefficients can be freely specified for 
So unless realization of t^^Cs) is more important than the realization 
t^iCs), P(s) should be chosen lower triangular. For this example 
P(s) is forced to be a lower triangular matrix. 

As before, assume that the transfer matrix of Eq, (A, 29) is still 
the desired transfer matrix but now the emphasis is on the poles of diago- 
nal terms. It is desired that the overall system poles be the poles of 
tfiCs) and Since P(s) is forced to be lower tri- 

angularjthe above is realized by the following choice of P(s) 


PCs) 


dll(s+2) Cs+3) 


d^iS 


*’ 22 ® 


+ h 


21 


0 


d22(s+^) 


C5.17) 


In the above d. . and h.. are as defined in Sten 4. 

11 ij ^ 

The desired K^’s are given by Eq. (4.30). But, as mentioned in 
Theorem 5.1, choice of the pole positions restricts the allowable free- 
dom of the K.'s as follows: 



12 /» 

1. only the elements in for which i ^ j can be chosen 
arbitrarily; 

2. for (i,j)th element in the transfer matrix, only additional 
error coefficients can be chosen arbitrarily for i > j * 


Thus, the specifiable part of the desired are given below 


ki - 


1 

1 

, k., = I 

* 

, k- = 


1 

1 1- 

2 

-7 *_ 

3 

-.39 

(5.18) 


where * indicates unspecif iable element. 

Step 7 : 

Since is nonsingular definite, this step is not necessary. 

Step 8: 

Since P(s) is forced to be a lower triangular matrix, Eq. (5.11) 
and (5.12) are used to solve for and K^. Also, since m=2, these 
equations need be solved only for i = 1,2 but in sequence. 

Solving Eq. C5.ll) for i 1 gives i = m-i+1 = 2 and hence 

^^ 22^1 "" ^^ 22^1 ^^ 22^1 

« (l.O)’^(l.O) = 1.0 

and solving Eq. j(5.12) for i = 1 gives I = m-i + 1 = 2. Hence 


0^12)1 = 


^^ 22^1 


{(ri2)i) 


Y {0} = 0 
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For i = 2, = m-i+l = 1, and hence from Eq . (5.11) 


Pll 

^21 


-1 


1 

0 


' 1 


Y 

1 

1 


10 


9 ^ 



k 



r. .. 


11 

12 


11 

1 

k2i 

Ir 

22 


^21 


Thus, and re completely known and are given as follows 


"l = 


1 d 
9 1 




1 0 
1 1 


(5.19) 


Step 9: ■ 


From Eq. (5.17), 


= P(s) 


s=o 


**21 ^**22 


If this matrix is compared with Eq. (5.19), the following is readily 
obtained 


21 


« 


‘^ll ^22 ^ 


9 
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Substituting in Eq. (5.17) completely specifies the diagonal terms of 
P(s) as follows 


P(s) = 


1/6 Cs+2)(s+3) 


0 


^21^^ ^22^ ^ l/4(s+4) 


Hence 


1 + 5/6 s + 1/6 s' 
9 + h22S + 


1 + 1/4 s 


p = ! 

1 0 


.83 0 ^ 


'.17 o' 



, P.. = 




1 ! 

9 1 

’ 2 

>2 

’ 3 

1 

_^21 


(5.20) 


Step 10: 

Once again, since P(s) is forced to be a lower triangular matrix, 

Eq. (5.13), (5.14) and (5.15) are used to determine the remaining 

unknown and unspecif iable components of P(s) and T(s), 

To find and K^, solve those equations for j=2. Eor i = 1, 

2 

I = m-i+1 = 2. But [p^] is completely known from Eq. (5.19) and, hence 
[k^J is found directly from Eq. (5.13). 
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Also, from Eq, (3.19) and (5.16) 



0 

0 


Hence , 




1 

1 


roi 


j 

0 
( 

i 

(.25) 






_.25_ 


For 1=2, £ = m-i+1 = 1 and hence from Eq. (5.14) 


(5.21) 


^^ 21^2 " ^^ 22^1 


(q2i)i + 


^ 22^2 


11 

^21 


^^2l\^^11^2 


Notice that (^22^2 l»een obtained in Eq. (5.21), (^ 2^)2 is the arbi- 
trary (desired) error coefficients specified in Eq. (5.18) and all 
other elements quantities on the right are known. Hence 

(P2j^>2 = (1)'^ (2 + [-7 .25] 

= 4.12 

Finally, the unknown part of is determined by solving Eq . (5.15) 
for i = 2. 


1 - 

,9. 


- (1)(.83) 
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For i = 2 , S/ — m-i+1 = i and hence 


<\l>2 “ 


= J {[I 0 ] 



"12^ 

^11 

11 

_^12_ 

~ .83 



>-12_ 

- (1) - 


- (q 


11 


^1 " ^^I2^2^^21^il 


.1 


= .83 - 1 = - .17 


The 


and 


matrices are known completely as given below 


\83 

0 ‘ 


-.17 

”1 

4.12 

.25 

«j 

’ 2 

■ 

0.7 

.25J 


(5.22) 


At last, proceed to find elements in and K^. To do that, first 
calculate as given by Eq. (3.19) 


Qa = R3 + K^Pa 


-.17 0 ' 


~ .83 o' 


[-.139 0 

• 

0.7 .25 _ 


4.12- .25^ 

! 

[1.613 .062^ 


(5.23) 


2 2 

For i = 1, Z ~ m-i+1 = 2 and since [p^] is knovm, [k^] is directly 
calclated by solving Eq. (5.13) for j = 3 
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[k^] 


1 


{[kj^] (?22^3 


[q^]> 



For i = 2, A = m-i+1 = 1 and the unspecified part of may be calcu- 
lated by solving Eq. (5.14) for j = 3 


^^ 21^3 ~ ^^ 22^1 t ^^ 21^2 ^^21 ^ 22^3 


11 


21 


^^21^1^^11^3 


L —-11 


= (1) ^ ^ 1.613 + [-.39 -.062] 


- (l)(.17)i 


= 1,613 - .39 - .5625 - .17 = .491 


Notice that the solution of the above equation was made possible by 
ensuring the inverse of (^22^1' in were so chosen as to 

force as obtained in Eq. (5.19) to be a lower definite matrix. 

Finally, after substitution of the results obtained so far, the 
only unknown part of is obtained by solving Eq. (5.15) for j = 3 
and i = 2. For 1=2, i - m-i+1 = 1 and hence 
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= Yitl 0] 


‘^11 

1 

Pll 

Pl2_ 

^17 ’ 

J 


_.491 

- (-.139) 




= .17 + ,139 = .31 


Thus, and are completely known as 


'’3 = 


’.17 

o] 



".31 

0 

.491 

. 


It 

-.39 

-.062 


(5.24) 


The design is almost compete, except that one has to find the 
pair {F,G} which corresponds to the components of P(s) obtained in 
Equations (5.19), (5.22), and (5,24). This can be trivially done as 
was done for the example (4.1). However, the first order of business 
is to check whether the transfer matrix so obtained is satisfactorily 
close to the one wanted. This is done as follows: 


T(s) = R(s)P(s) 


-1 


s+1 

0 ^ 


,2(s+5) 

ij 

L 

fs+1 

0 ' 

n 

[ 2 ( 3 + 5 ) 

1 



. 17 (s + 5s + 6) 


.491s + 4.12s + 9 


,25(s+4) 


0 

.25(s+4) 

0 


“1 


-(.491s“ + 4.12s + 9) .17(s + 5s + 6) ! 


/A 



131 


where A = 1/24 (s+2) (s+3) (s+4) . 


Hence 


T(s) 


6(s+l) 

(s+2)(s+3) 

24(.009s^ + .38s + 1) 
(s+2) (s+3) (s+4) 


0 

4 

s+4 


6(s+l) 
(s+2) (s+3) 


.216s^ + 9.02S+24 4 

( s+2 ) ( s+3 ) ( s+4 ) s+4 


(5.25) 


Once again the transfer matrix obtained above could be compared 
with the one desired as given in Eq, (4.29), and repeated here for 
convenience . 


T(s) 


6(s+l) 
(s+2) ( s+3) 

IQ 

(s+2) (s+5) 


0. 6(3+1) 
(s+2) (s+3) 

4 

s+4 


(5.26) 


Since the emphasis is on the poles of the diagonal terms, one 
can see they are exactly at the desired place. However, this is possible 
only by the sacrifice of the term (5.26). The situations 

in which the off diagonal terms come out to be zero is a special case, 
as is discussed in the next section. The term transfer ma- 

trix of Eq. (5.25) has the same generalized error coefficients as speci- 
fied by the corresponding term in desired transfer matrix in Eq . (5.26). 
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However, the pole positions of t^^^Cs) are not the same as desired. It 
is up to the designer to test this transfer function to see if its 
transient response is satisfactory and within his tolerance limits. 

The step responses of the two are compared in fig. 5.1, 

If the transfer matrix obtained above is satisfactory, the 
pair {F,G} necessary to realize it are easily found by proceeding on 
exactly the same lines as was done for example 4.1. Here, only the 
results are given for the sake of completeness. 



1.5 


- 2.21 


-.25 


-2.19 


0 

-2 


G = 


6 


[- 11.8 


o' 

4 


(5.27) 


In the next section some special case are discussed. 

5.6 Special Cases 

5.6.1 The Case of Triangular Decoupling 

It is seen in Example 5.1 that by fixing the poles in advance, 
the resulting transfer matrix came out to be lov/er triangular, even 
though it was not intended so. This may have certain advantages if it 
could be predicted in advance. The following lemma gives the suffi- 
ciency condition for triangular decoupling . 



IT STE 
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Lerana 5.1 If R(s) is upper (lower) triangular, then T(s) 
is upper (lower) triangular decoupled by choosing P(s) 
upper (lower) triangular, provided = R(s) is 

nonsingular. 

-1 

Proof: If P(s) is upper (lower) triangular, then P (s) 

is also upper (lower) triangular, and hence T(s) = R(s) 

P ^(s) is also upper (lower) triangular. This together 
with Theorem 5.1 proves the lemma because R(s) triangular 
and R^ nonsingular implies R^ is absolutely nonsingular 
[Definition 5.1 in Section 5.2], 

The main advantage gained by above is that instead of having no control 
over the one side of the diagonal transfer functions (Theorem 5.2), 
they can be forced to be zero. 

Since R(s) is completely known in advance before any actual 
design is attempted, it could be easily found whether R(s) is a tri- 
angular matrix or can be made triangular by simple changes in the names 
of inputs and outputs. Particularly if the open loop system is coupled 
in one direction only, then a change in the names of inputs and out- 
puts would make-R(s) a triangular matrix. Of course, a change in the 
output matrix can cause a significant change in R(s) . If such a change 
is not intolerable, then it might be worth forcing R(s) triangular. 

Thus, if R(s) is triangular, then the designer knows in advance 
that by using the pole fixing method, a triangular transfer matrix design 
can be obtained. This has distinct advantage over conventional decoupled 
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design where only diagonal transfer functions are nonzero and all other 
transfer functions are forced to zero, 

5.6.2 The Minimum Constraint Case 

The case in which there is a minimal constraint on the general- 
ized error coefficient is stated through the following corollary. 

Corollary 5.1: If all a’s are equal, i.e. , “ ^2 ” 

= o = then for each transfer function in the transfer 
m 

matrix T(s) , a + 1 arbitrary generalized error coefficients 
can be realized if : 

1. is chosen nonsingular. (This is a necessary condition 
for the system to be stable and to have nonsingular transfer 
matrix as stated in Theorem 3.1.) 

2. K , - is chosen so that P is nonsingular. 

Proof: The proof follows from Theorem 4.1. In addition, note 

that for any choice of P , , , K , > is uniquely determined via 
Eq. (3.19), since and, hence, P^ are nonsingular. This 
guarantees that there are many sets of which result in 

nonsingular 

This case of equal o*s is mentioned here for its academic 
value only. In practice, the o's, even though non-unique are feed- 
back invariants. Thus, if n is not an integral multiple of m, it is 
impossible to make them equal without adding dynamics to the system. 
The case in which dynamics are added to the system is not considered 
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In this study and in fact, has been considered by others, the most 
recent of which is the work of Wolowich [1973] . 

5 . 7 Summary 

In this chapter a new design method called the pole fixing 
method is used for the synthesis of multivariable systems. In many 
cases the pole fixing method can be successfully used for the design 
of multivariable systems where the conventional design methods or 
the error coefficient method of the last chapter fails. In this method 
all the poles of the multivariable system are pre-assigned and the 
zeros are adjusted to produce satisfactory steady state behavior. The 
only disadvantage of the pole fixing method is that it puts new con- 
straints on the feedback invariant part of the system, R(s), But it is 
seen that in many cases, simple changes in the output matrix satisfy the 
artificial constraints introduced by the pole fixing design method. 



CHAPTER 6 


SUMMARY AND CONCLUSION 

This chapter suimnarizes the results of this study and pertinent 
results concerning linear state variable feedback (l.s.v.f.) invariants. 
Some areas of further research in connection with this study are also 
indicated. 

6.1 Summary 

Linear state variable feedback has been used in the design of 
multivariable systems. Multivariable systems differ from the scalar 
case in that l.s.v.f. affects not only the pole positions, but also 
directly affects the zeros because of the coupling between the sub- 
systems. 

The l.s.v.f. invariants of multivariable systems were not known 
completely until recently. Wolowich [1972] introduced a complete set 
of l.s.v.f. invariants for multivariable systems for the first time. 

He utilized the structure of multivariable systems to derive l.s.v.f, 
invariants and, in turn, used l.s.v.f. invariants to derive the neces- 
sary and sufficient condition for exact model matching. If l.s.v.f, 
can be found, such that the closed loop response matches a pre- 
specified (model) response, then the plant (open loop system) is 
said to match with the model. 
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The necessary and sufficient conditions for exact model matching 
did not come as a relief to the designer. The shortcoming is that if 
the necessary and sufficient conditions are not met, there is no way to 
complete the design. One does not know what changes need to be made 
in the design specification to satisfy the necessary and sufficient con- 
ditions. Thus, the problem of designing a multivariable system using 
l.s.v.f. where noninteraction is not a design criteria remained largely 
unanswered. An exception is a very special case where the dynamics of 
each element in any given row of the transfer matrix are the same except 
for the gain [ferg, 1971], 

In this study two complete design procedures are developed for 
the first time for the design of multivariable systems incorporating 
cross coupling using l.s.v.f. alone. Noninteraction is treated as a 
special case of interaction. The main design equation is developed in 
Chapter 3. The first design procedure using error coefficients alone 
as the design criteria is worked out in Chapter 4, This is followed 
by a simple example to illustrate the method. 

The generalized error coefficients represent the system input- 
output relationship only after the system transients have decayed to 
insignificant amplitudes. Here lies the shortcoming of the error co- 
efficient design method, since it does not take into account the 
transient behavior of the system. To alleviate this problem the pole 
fixing method is introduced. As the name suggests, all the poles of 
the system are fixed in advance ensuring system stability and fast 
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decaying transients. If the system satisfies certain sufficiency 
conditions, then in addition to fixing all the poles at desired places, 
some indirect control can be exercised over zeros. In addition, those 
systems which cannot be decoupled can be designed using the pole fixing 
method. Obviously, this method of design does not produce any uncontrol- 
lable poles as does the noninteracting design. 

The sufficiency conditions are variant under the change in sys- 
tem output matrix. Thus, in many cases a simple operation like changing 
the name of inputs or outputs could be enough to satisfy the sufficiency 
conditions. The application of the pole fixing method is illustrated by 
working a simple example. 

6 . 2 Fur ther Res ear ch 

Although this study provides two design methods for the design 
of multivariable systems. The following related topics merit further 
research: 

1. Establish additional definite relationships between the speci- 
fiable error coefficients, the poles of the overall system, and 
the component matrices of the matrix P(s). The first of the 
two design methods for the design of multivariable systems is 
based upon the assumption that the system input-output behavior 
can be satisfactorily described in terms of finite error co- 
efficients. In many cases this may not be true. In particular, 
when a transfer function is expanded in terms of error coeffi- 
cients, much of the information about the transient response is 
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lost because only a relatively small number of error coeffi- 
cients are considered. The transfer function cannot always be 
uniquely reconstructed only from the specification of these 
error coefficients. 

2. Investigate the relation between the components of a polynomial 
matrix and its determinant. This is useful because poles of 
the overall system are given by the zeros of determinant 

P(s). The components of P(s) are determined one by one. The 
relation between the components of P(s) and its determinant 
. can then be used to determine whether the particular component 
of P(s) is suitable enough for the design to proceed any fur- 
ther. If not, the error coefficients could be suitably 
changed before the design is complete. For example, the 
constant term in the determinant of P(s) is given by the 
determinant of P^, the first component of P(s). 

3. Partially or completely fix the poles of the systems without 
forcing P(s) a triangular matrix. By forcing P(s) a triangular 
matrix the designer could of course fix all the poles of the 
system in advance. But this has the disadvantage that control 
is lost, over the zeros of one side of diagonal terms in T(s). 
Also, fixing all poles adds additional constraints which must 
be satisfied. Thus, forcing P(s) to be triangular might be 
over restrictive. Two obvious ways to proceed are to either 
partially fix the poles or equivalently to find a way in which 
P(s) is not forced to be a triangular matrix. 
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4. Find the necessary condition for the design method to be 
applicable. Theorem 4.1 and 5.1 provide only sufficient 
conditions to complete the design. Upon investigation it 
might turn out that some of the sufficiency conditions are 
also necessary conditions. The advantage of knowing the 
necessary conditions is that the designer knows that he 
is not being too conservative in designing the system. 

( 

It is seen that even though linear state variable feedback 
(l.s.v.f.) can be successfully used for the. design of multivariable 
systems, all of its implications in the design of multivariable feed- 
back have not been resolved. It is expected that efforts of more than 
just a few will be required before those working in the area feel that 
the problem is solved. 
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